Small time Edgeworth-type expansions for weakly convergent 
nonhomogeneous Markov chains.* 

Valentin KONAKOV 
Central Economics Mathematical Institute, Academy of Sciences 
Nahimovskii av. 47, 117418 Moscow, Russia 
E mail: kv24@mail.ru 
Enno MAMMEN 
Department of Economics, University of Mannheim 
L 7, 3-5, 68131 Mannheim, Germany 
E mail: emammen@rumms.uni-mannheim.de 

February 1, 2008 



Abstract 
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1 Introduction. 



Recently, there was some activity on Edgeworth-type expansions for dependent data. In most approaches 
higher order expansions have been derived by appHcation of classical Edgeworth expansions for inde- 
pendent data. The approaches differ in their main idea how the dependence structure can be reduced 
to the case of independent data. For sums of independent random variables and for functionals of such 
sums the theory of Edgeworth expansions is classical and well understood in a very general setting (see 
Bhattcharya and Rao (1976) and Gotze (1989)). For models with dependent variables three approaches 
have been developed where the expansion is derived from models with sums of independent random vari- 
ables. In the first approach mixing properties are used to approximate the Markov chain by a sum of 
independent random variables and it is shown that their Edgeworth expansion carries over to the Markov 
chain up to a certain accuracy. The mixing approach was first used by Gotze and Hipp (1983) and it 
was further applied to continuous time processes in Kusuoka and Yoshida (2000) and Yosliida (2004). 
Under appropriate conditions Markov chains can be splitted at regeneration times into a sequence of i.i.d. 
variables. This fact has been used in Bolthausen (1980, 1982) to get Berry-Esseen bounds for Markov 
chains. For the statement of Edgeworth expansions the regenerative method has been used in Mahnovskii 
(1987), Jensen (1989), Bertail and Clemencon (2004) and Fukasawa (2006a). The higher order Edgeworth 
expansions have been used to show higher order accuracy of different bootstrap schemes, see Mykland 
(1992), Bertail and Clemencon (2006) and Fukasawa (2006b). 

Both approaches, the mixing method and the regenerative method only have been used for Markov 
chains with a Gaussian Hmit. In this paper we study Markov chains that converge weakly to a diffusion 
Hmit. For the treatment of this case we make use of the parametrix method. In this approach the 
transition density is represented as a nested sum of functionals of densities of sums of independent 
variables. Plugging Edgeworth expansions into this representation will result in an expansion for the 
transition density. Thus as in the mixing method and in the regenerative method the expansion is 
reduced to models with sums of independent random variables. 

The parametrix method permits to obtain tractable representations of transition densities of diffusions 
and of Markov chains. For diffusions the parametrix expansion is based on Gaussian densities, see Lemma 
1 below, and standard references for the parametrix method are the books of Friedman (1964) and 
Ladyzenskaja, Solonnikov and Ural'ceva (1968) on parabolic PDE [see also McKean and Singer (1967)]. 
For a short exposition of the parametrix method, see Section 3 and Konakov and Mammen (2000). 
Similar representations hold for discrete time Markov chains Xkji, see Lemma 3 below. The parametrix 
method for Markov chains was developed in Konakov and Mammen (2000) and it is exposed in Section 
3.2. In Konakov and Mammen (2002) the approach was used to state Edgeworth-type expansions for 
Euler schemes for stochastic differential equations. Related treatments of Euler schemes can be found in 
Bally and Talay (1996 a,b), Protter and Talay (1997), Jacod and Protter (1998), Jacod (2004), Jacod, 
Kurtz, Meleard and Protter (2005) and Guyon (2006). 

In this paper we study triangular arrays of Markov chains Xk^h [k > 0) that converge weakly to 
a diffusion process Yg (s > 0) for n — > oo. We consider the Markov chains for the time interval 
(0 < A; < n). The corresponding time interval of the diffusion is (0 < s < T). The term h = T/n denotes 
the discretization step. We allow that T depends on n. In particular, we consider the case that T — > 
for n — !■ oo. Furthermore, we allow nonliomogeneous diffusion limits. 

Weak convergence of the distribution of scaled discrete time Markov processes to diffusions has been 
extensively studied in the literature ( see Skorohod (1965) and Stroock and Varadhan (1979)). Local limit 
theorems for Markov chains were given in Konakov and Molchanov (1984) and Konakov and Mammen 
(2000, 2002). In Konakov and Mammen (2000) it was shown that the transition density of a Markov 
chain converges with rate 0(n~^/^) to the transition density in the diffusion model. For the proof there 
an analytical approach was chosen that made essential use of the parametrix method. 

The main result of this paper will give Edgeworth type expansions for the transition densities of the 
Markov chains X^^h (0 < k < n). The first order term of the expansion is the transition density of the 
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diffusion process (0 < s < T). The order of the expansion is o{hr^^^) with 6 > 0. Related results 
were shown in Konakov and Mammen (2005) . The work of this paper generalizes the results in Konakov 
and Mammen (2005) in two directions. The time horizon T is allowed to converge to and also cases 
are treated with nonhomogeneous diffusion limit. Small time asymptotics is done for two reasons. First 
of all it allows approximations for the joint density of values of the Markov chain at an increasing grid 
of points. Secondly, it is motivated by statistical applications. In statistics, diffusion models are used 
as an approximation to the truth. They can be motivated by a high frequency Markov chain that is 
assumed to run in the background on a very fine time grid and is only observed on a coarser grid. If 
the number of time steps between two observed values of the process converges to infinity this allows 
diffusion approximations (under appropriate conditions). This asymptotics reflects a set up occurring 
in the high frequency statistical analysis for financial data where diffusion approximations are used only 
for coarser time scales. For the finest scale discrete pattern in the price processes become transparent 
and do not allow diffusion approximations. The statistical implications of our result will be discussed 
elsewhere. The mathematical treatment of nonhomogeneous diffusion limits with time horizon T going to 
zero contributes some additional qualitatively new problems. In this case some additional terms appear 
that explode for T ^ and for this reason these terms need a qualitatively different treatment as in the 
case with fixed T. The nonhomogeneity adds an additional term in the Edgeworth expansion. See also 
below for more details. 

The paper is organized as follows. In the next section we will present our model for the Markov chain 
and state our main result that gives an Edgeworth-type expansion for Markov chains. Connections with 
previously known results are also discussed in Section 2. In Section 3.1 we will give a short introduction 
into the parametrix method for diffusions. In Section 3.2 we will recall the parametrix approach developed 
in Konakov and Mammen (2000) for Markov chains. Technical discussions, auxiliary results and proofs 
are given in Sections 4 and 5. 

2 The main result: an Edgeworth-type expansion for Markov 
chains converging to diffusions. 

We consider a family of Markov processes in R'^ that have the following form 

Xk+i,h = Xk,h + m {kh, Xk,h) h + y/h^k+i,h, ^o,/» = a; e R'', fc = 0, n - 1. (1) 

The innovation sequence (Ci^/j) ^ is assumed to satisfy the Markov assumption: the conditional 
distribution of (,k+i h given the past Xk^^ = Xk, ...,Xo,/i = xq depends only on the last value Xk^h = Xk 
and has a conditional density q {kh, Xk -r). The conditional covariance matrix corresponding to this density 
is denoted by a{kh,Xk) and the conditional v — th cumulant by Xv{^h,Xk)- The transition densities of 
{Xi^h)i^\ n ^-re denoted by ph (0, kh,x, •). The time horizon T = T{n) < 1 is allowed to depend on n 
and h = T/n. is the discretization step. 
We make the following assumptions. 

(Al) It holds that /^<j yq{t,x,y) dy = for < t < 1, x€ W^. 

(A2) There exist positive constants cr* and a* such that the covariance matrix a {t, x) 
satisfies 

(T* < e^a [t, x)e<a* 
for all 116111 = 1 and i e [0, 1] and x e K'^. 



= /iRd yy^Q {t, X, y) dy 
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(A3) There exist a positive integer S' and a real nonnegative function tp (y) ,y gW^ satisfying supj^ggd ip {y) < 
00 and /gd V {y) dy <oo with S = {S' + 2)d + 4 such that 

\Dlq{t,x,y)\<^{y), [0,1], x,y eW^ |H =0,1,2,3,4 

and 

\Dlq {t, x,y)\<i;{y), te [0, 1], x,yGR'^ = 0, 1, 2. 

Moreover, for all x,y £ R'^, h > 0,0 < t,t + jh < l,j > jo, with a bound jo that does not depend 
on x, t, 

D^q^J) {t, x,y)\< Cj-"'^^ (j-"^y) , = 0, 1, 2, 3 

for a constant C < oo. Here q^^^ {t, x, y) denotes the j-fold convolution of q for fixed a: as a function 
of y: 

q'^\t,x,y) = j q'^^~'^\t,x,u}q{t+{j -l)h,x,y-u)du, 

q^'^\t,x,y) = q{t,x,y). 

Note that the last condition is motivated by (A2) and the classical local limit theorem. Note also that 

for 1 < j < jo 

" \y\f q^^\t,x,y)dy<C{jo,S). 
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(Bl) The functions m {t, x) and a {t, x) and their first and second derivatives w.r.t. t and their derivatives 
up to the order six w.r.t. x are continuous and bounded uniformly in t and x. All these functions 
are Lipschitz continuous with respect to x with a Lipschitz constant that does not depend on t. The 
functions X}At:^)j W\ = 3,4, are Lipschitz continuous with respect to t with a Lipschitz constant 
that does not depend on x. A sufficient condition for this is the following inequality 

/ (1 + jj^jj*) \q {t, x,z)-q {t', x,z)\dz <C\t-t'\,0<t,t' <1, 

with a constant that does not depend on a; G W^. Furthemore, D^a{t,x) exist for < 6 and are 
Holder continuous w.r.t. x with a positive exponent and a constant that does not depend on t. 

(B2) There exists x < ^ such that lim inf„^oo 7'(n)n'^ > 0. 

The Markov chain X^^^, see (1), is an approximation to the following stochastic diflterential equation 
in M'^ : 

dYs = m {s, Y,) ds + A {s, Y^) dWs, Yq = x e M!^, s e [0, T], 

where (W^s)s>o standard Wiener process and A is a symmetric positive definite dx d matrix such 

that A{s,y) A(s,y)'^ = a{s,y) . The conditional density of Yt, given Iq = a; is denoted by p{0,t,x,-). 
We will use the following diflferential operators L and L : 

Lf{s, t,x,y) = ^Yl '^ij ^) ^ '^dx-dx-^^ +^mi{s,x) ^"^^ g^'."^' , 

i,j=l * i=l * 

Lf{s,t,x,y) = - 2^ aij{s,y) ^^^^ +}^mi{s,y) — . (2) 

i,_7 = l * i=l * 
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To formulate our main result we need also the following operators 

, daij{s,x) d^f{s,t,x,y) ^ ^ dmi{s,x) df {s,t,x,y) 



~, 1 daij{s,y) d"^/ {s,t,x,y) dmi{s,y) df {s,t,x,y) 
Lf{s,t,v,z) = -^ —g-^ ___+^^_ _ . 



dxidxj 



i=l 



ds 



dxi 



(3) 



and the convolution type binary operation ® : 

f I® g{s,t,x,y) = I du f {s,u,x,z)g{u,t,z,y)dz. 

Konakov and Mammen (2000) obtained a nonuniform rate of convergence for the difference ph (0, T,x,-)- 
p (0, T, X, •) as n — > oo in the case T x 1. Edgeworth type expansions for the case T >c 1 and homogenous 
diffusions were obtained in Konakov and Mammen (2005). The goal of the present paper is to obtain 
an Edgeworth type expansion for nonhomegenous case which remains vahd for the both cases T x 1 or 
T = o(l). The following theorem contains our main result. It gives Edgeworth type expansions for Ph- 
For the statement of the theorem we introduce the following differential operators 



J^i[f]{s,t,x,y) = 



kl=3 



^^All^D^J{s,t,x,y), 
J^2[f]{s,t,x,y)= J2 ^^^^D:f{s,t,x,y). 

1^1=4 

Furthermore, we introduce two terms corresponding to the classical Edgeworth expansion (see Bhat- 
tacharya and Rao (1976)) 



Tri{s,t,x,y) 
T^2{s,t,x,y) 



{t-s)J2 ^^^'^^Dlp{s,t,x,y), 
(t-s) E ^"^"'l'^^ Dlp{s,t,x,y) 



(4) 



^Ut-sflY^^^^^D':^) pis,t,^,y), 

.kl=3 



(5) 



where 



xAs,t,y) 



t - s 



xAii,y)du 



and Xu{tj^) is the u — th cumulant of the density of the innovations q{t,x,-). The gaussian transi- 
tion densities p{s,t,x,y) are defined in (6). Note, that in the homogenous case Xvi'^iV) = Xuiv) ^^'^ 
X^(s,f, y) = Xv{y)j where Xviv) is the v — th cumulant of the density q{y, •). 

Theorem 1. Assume (Al)-(A3), (B1)-(B2). Then there exists a constant 5 > Q such that the following 
expansion holds: 

s'\ 



sup T''/^ I _^ 



y ^ X 



Vt 



Ph{0,T,x,y)-p{0,T,x,y) 
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-/i^/^TTi (0, T, X, y) - hTT2{0, T, X, y) 



= 0{h'+% 



where S'is defined in Assumption (A3) and where 

7ri{0,T,x,y) = {p ^ J^i\p]){0,T,x,y), 

7T2{0,T,x,y) = {p^J^2\p]){0,T,x,y)+p^J^i\p^J^i\p\]iO,T,x,y) 

+ ^p ® {Ll - L2)p(o, T, x, y) - ® (L' - L')p(0, T, x, y). 

Here p{s,t,x,y) is the transition density of the limiting diffusion Ysand the operator is defined as L, 
but with the coefficients "frozen" at the point x. The norm ||-|| is the usual Euclidean norm. 

Remark 1. The terms of the Edgeworth expansion have subgaussian tails and are of order n~^/^ or 
n^^, respectively: 



h^/^TTi{0,T,x,y) < Cin-i/^T-'^/^exp 
|/i7r2(0,T,a;,t/)| < Cin-^T-'^/^ exp 



-a 



-Co 



y-x 



Vt 

y-x 



VT 



with some positive constants Ci and C2. 



Remark 2. If the innovation density q{t,x,-) and the conditional mean m{t,x) do not depend on 
X then we are in the classical case of independent non identically distributed random vectors. We now 
show that then the Edgeworth expansion of Theorem 1 coincides with the first two terms of the classical 
Edgeworth expansion ft,^/^7ri(0, T, x, y) + hTT2{0, T, x, y). Note first that in this case = L,L' = L' and 
p{s,t,x,y) = p{s,t,x,y) where p is defined in (6) with a{s,t,y) — cr{s,t) = a{u)du and m{s,t,y) = 
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m(s, t) = m{u)du. This gives 



7ri(0,T,a;,2/) = / ds I p{0,s,x,v) ^ ^^^^^Dlp{s,T,v,y)dv 
^^"^^dsDy I p(0,s,x,v)p{s,T,v,y)dv 



E 

kl=3 



J2 %.(0,T)£>;p(0,r,ar,y) 



1^1=3 

T 



^ -X. (0, T)£>^p(0, T, y) = 5?i (0, T, x, y), 



kl=3 



p(^J^i\p\{s,T,z,y)= f du f p{s,u,z,w) ^ ^^^^^D!^p{u,T,w,y)dw 

•''^ |i/|=3 



dwi)^p(s,T,0,y) = (T-s) ^ 



kl=3 



kl=3 



E 



(^-^) E 



,k'l=3 

D-+-''p{s,T,z,y), 



IH=3,k'l=3 



+ f ds j p{0,s,x,z){T-s) ^^''^^D'^z^''Pi''T,z,y)dz 

•'° •' |H=3>'|=3 ^' ^' 



kl=4 



= ^E 



J/! 



|,.|=4 

For = z^' we have 
Jo 

For y ^ v' we get 



|i,|=3,|^'|=3 ■ Jo \Js J 



T rT 



X^'{u)du\ ds = - 



xAs)xAu)dsdu = — x,(0,T)x,(0,T). 







X^(s) Xv'{u)du ]ds+ I x,.'(s) I / Xu{u)du \ ds 

IQ Js 



^ L Is ^'^'^ ^ ^^^^^ ^"''^ '^^^^ 

= \j^ / [Xi/(s)x,.'(w) + Xu'{s)xA'U')\ dsdu 
= ^X.(0, T)x,, (0, T) + ^x.' (0, T)x,(0, T). 
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From these equations we obtain 



P «) [p] (0, T, a;, y) + p [p ® jTi [p] (0, T, a;, y) 




= Tj2^^D:mT,x,y) + ^lY. 



Dl) mT,x,y) 




This shows the claim that we get for this case tiie first two terms of the classical Edgewortli expansion. 

Remark 3. If xAt,x) = for \u\ = 3 and for t € [0,T] x R"^ then it holds that J^i = 0. The Theorem 
1 holds with 



If in addition XuitiX) = for = 4 then the first four moments of the innovations coincide with the 
first four moments of a normal distribution with zero mean and covariance matrix a{t, x). In this case we 
have = and we have 



and the first two terms of the Edgeworth expansion do not depend on the innovation density. In particular, 
it holds that Xy(t, .t) = for = 3,4 for Markov chains that are defined by Euler approximations to 
diffusions. Thus, an Edgeworth expansion for the Euler scheme holds with the same tti and 7r2 as just 
defined. For the homogenous case we have that L' = L' = and we obtain for the Euler scheme in this 
case 



This result for T = [0, 1] under Hormander's condition on a diffusion matrix was obtained by Bally and 
Talay (1996). 

Remark 4. We now shortly discuss an application of our result to statistics. Assume that one observes 
a Markov process Xi^h, Xnk.h at time points k, 2k, nk. That means we assume that a high frequency 
Markov chain runs in the background on a very fine time grid but that it is only observed on a coarser grid. 
This asymptotics reflects a set up occurring in the high frequency statistical analysis for financial data 
where diffusion approximations are used only for coarser time scales. For the finest scale discrete pattern 
in the price processes become transparent that could not be modeled by diffusions. The joint distribution 
of the observed values of the Markov process is denoted by Ph- We assume that this joint distribution 
can be approximated by the distribution of (Yi, Yn) where Yi, Yn are the values of a diffusion on the 
equidistant grid fc/i, 2kh, nkh. The joint distribution of (Yi, Y„) is denoted by Qh- According to our 
theorem the one-dimensional marginal distributions of Ph, can be approximated by the one-dimensional 
marginal distributions of Qh- Under appropriate conditions the Li-norm of this difference is of order 
^-i/2_ 'jjjjg iniplies that the Li-norm of the difference between the joint distributions Ph and Qh is of 
order nk~^^^ . That means the diffusion approximation is only accurate if /c ^ n^, i.e. only if the grid of 



7ri(0,T,a;,y) = 0, 

772(0, T, X, y) = {p® T2\p]m T, X, 2/) + ® (Ll - L'')p{0, T, x, y)-^p® (L' - L')p{0, T, x, y). 



7ri(0,T, x,y) 
7r2{0,T,x,y) 



0, 

^p g) {Ll - L')p{Q, T, x, y)-^p® {L' - L')p{0, T, x, y) 



ni{0,T,x,y) = 0, 

n2{0,T,x,y) = ^p® {L^- L'')p{0,T,x,y). 
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observed points is verj' coarse in comparison to the grid on which the Markov process lives. Only in this 
case it can be guarantied that a statistical inference that is based on the diffusion model is accurate. Or 
put it in another way, data that come from the Markov model could not be asymptotically statistically 
distinguished from diffusion observations. Our results help to analyze what may go wrong if fc does 
not hold. The (signed) transition densities j? + /i^/^7ri + /i7r2 given in the statement of Theorem 1 define a 
joined (signed) measure Rh- According to Theorem 1, the marginal distributions of Rh approximate the 
onc-dimcnsional marginal distributions of Ph with order o(fc^^^*). One may conjecture that under some 
regularity assumptions the exact order is k^'^/'^ . This imphes that \\Ph — Rh\\i is of order nk~^^^. Thus, 
this approximation is appropriate as long as fc » n^/^. This is a much more acceptable assumption. 
Now, one can check which statistical procedures behave differently under the models Qh and Rh- These 
procedures may lead to erroneous conclusions for the Markov data. 

3 The parametrix method. 

3.1 The parametrix method for diffusions. 

We now give a short overview on the parametrix method for diffusions. For any s €: [0,T], x,y € M.'^ we 
consider the following family of "frozen" diffusion processes 



Let {s,t,x, •) be the conditional density of Yj, given Yg = x. In the sequel for any z we will denote 
p{s,t,x, z) = (s,t,x,z) , where the variable z acts here twice: as the argument of the density and as 
defining quantity of the process Yf. 

The transition densities p can be computed explicitly 



dYt = m {t, y)dt + A {t, y) dWt, Y, = x, s 



<t<T. 



p{s,t,x,y) = (27r) '^^^ (deta (s, i,t/)) ^ 



X exp 



"2 (y~2;-TO(s,i,y)) a {s,t,y) {y - x - m {s,t,y)) 



(6) 



where 




m {u, y) du. 



Note that the following differential operators L and L correspond to the infinitesimal operators of Y or 
of the frozen process Y , respectively, i.e. 



Lf{s, t, X, y) = lim h-\E[f{s, t, Y{s + h),y)\ Y{s) =x]- f{s, t, x, y)}, 

n— >0 

Lf{s, t, X, y) = lim h-\E[f{s, t, F (s + h),y)\ Y{s) = x] - f{s, t, x, y)}. 

n— >0 



We put 



H={L- L)p. 



Then 




-o-y {s,y)) 



'p {s,t,x,y) 
dxidxj 
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In the following lemmas the k- fold convolution of H is denoted by H^''^ . The following results have been 
proved in Konakov and Mammen (2000). 



Lemma 1. Let < s < t < T. It holds 

OO 

p{s,t,x,y) = H^''\s,t,x,y). 



r=0 



Lemma 2. Let < s < t < T . There are constants C and Ci such that 

\H{s, t, x,y)\< Cip~^(t)c,p{y - 

and 

p® g^'-) {s, t, x,y)\< ^cAv - 

where =t — s, (f)c,p{u) = p~'^(l)c{u/ p)and 

eM-C\\ut) 



(j)c{u) = 



JoM-C\\vfdv) 



3.2 The parametrix method for Markov chains. 

We now give a short overview on the parametrix method for Markov chains. This theory was developed 
in Konakov and Mammen (2000). For any < jh < T, x,y G we consider an additional family of 
"frozen" Markov chains defined for jh < ih < T as 

Xi+i,h = Xi^h + m {ih, y)h + Vft^^+i,;,, Xj^h = x j < i < n, (7) 

where Cj+i,^' ■■■'^n,h innovation sequence such that the conditional density of ^j+i^/j given the past 
Xi^h = Xi, Xo^h = Xo equals to q{ih,y,). Let us introduce the infinitesimal operators corresponding 
to Markov chains (1) and (7) respectively, 

Lhf {jh, kh, x,y) = h~^ (^J Ph {jh, {j + 1) h, x, z) f {{j + 1) h, kh, z, y)dz- f {{j + 1) h, kh, x, y)j 

and 

Lhf {jh, kh, X, y) = h~^ (^j p^ {jh, {j + 1) h, x, z) f {{j + 1) h, kh, z, y)dz- f {{j + 1) h, kh, x, y)j , 

where {jh,j'h,x,-) denotes the conditional density of Xj/^h given Xj^h = x. Similarly as above, for 
brevity for any z wc write ph {jh, j'h,x, z) = pf^ {jh, j'h. x, z) . where the variable z acts here twice: as 
the argument of the density and as defining quantity of the process Xi^h- For technical convenience the 
terms / {{j + 1) h, kh, z, y) on the right hand side of Lhf and Lhf appear instead of / {jh, kh, z, y) . 
In analogy with the definition of H we put, for k > j. 



Hh {jh, kh, X, y) = {Lh - Lh^ Ph {jh, kh, x, y) . 
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We also shall use the convolution type binary operation iSih which is a discrete version of 0: 
g'^hf {jh, kh,x,y) = y2^ 9 {jh, ih, x, z) f {ih, kh, z, y) dz, 

where < j < k < n. We write g = g and g (g)^ Hj^^ — (^g (g)^ H/^^^^^ Hh for r > 1. For the 

higher order convolutions we use the convention Yl\^j = for i < j. One can show the following analog 
of the "parametrix " expansion for ph [see Konakov and Mammen (2000)]. 

Lemma 3. Let < jh < kh < T. It holds 

k-j 



Ph{jh,kh,x,y) = (g/i Hf^'ijh, kh,x,y), 



r=0 

where 



Ph{jh,jh,x,y) =ph{kh,kh,x,y) = 5{y - x) 
and 6 is the Dirac delta symbol. 

4 Some technical tools. 

4.1 Plugged in Edgeworth expansions for independent observations. 

In this Section we will develop some tools that are helpful for the comparison of the expansion of p (see 
Lemma 1) and the expansion oi Ph ( see Lemma 3). These expansions are simple expressions in p or 
Ph, respectively. Recall that p is a Gaussian density, see (6), and that ph is the density of a sum of 
independent variables. The densities p and ph can be compared by application of the classical Edgeworth 
expansions. This is done in Lemma 5 and this is the essential step for the comparison of the expansions 
of p and Ph- Lemmas 4 and 7 contain technical tools that will be used below. Lemma 7 contains bounds 
on derivatives of p;i that will be used at several places in the proof of Theorem 1. Its proof makes use of 
Lemma 6 that is a generalisation of a result in Konakov and Molchanov (1984) (Lemma 4 on page 68). 
Lemma 5 is a higher order extension of the results from Section 3.3 in Konakov and Mammen (2000). 

For the formulation of the lemmas we need some additional notations. Suppose that X e M'* is a 
random vector having a density q{x),x gR'', EX = 0,Cov{X,X) = S, where S be a positively definite 
d X d matrix . Denote A = ||ay || = and let Xv{^) be a cumulant of the order v = {vi, Vd) of 

a random vector Z G W^, <j){x) denotes a function in M!^ such that D^(f){x) exist and are continuous for 

=4, and A"^ = ||a*J'|| = T,^/^. 

Lemma 4. The following relation holds for s = 3 and for s = 4 

2^ lA 2^ ,A 



where z = Ax. 
Denote 



\u\=s |l/|=S 



k-1 k-1 



Hkiv) =h^m{ih,y),Vj,k{y) = h^a{ih,y). (8) 
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Proof of Lemma 4. For \u\ = 3,iy = (vi, Vd), each cumulant x,^(v4X) is a linoar combination of 
X^(^) with \n\ = 3 and with coefficients depending only on a^. It follows from the following relation 

Xvi^X) = iJ.„iAX) = J {aiixi + ... + aidXdY^ x ••• x (aiia;i + ... + auXdY'^q{'>^)dy.. 

Analogously, from the usual differentiation rule of a composite function and from the relation (f){z) = 
(f){Ax) , X = A~^z, it follows that D'^<f){z) = D'^(j){Ax) is a linear combination of D^(j){Ax) with coefficients 
depending only on a*^ . As a result of such substitutions we obtain that 



^ u\ 3! ^ 



|H=3 



E 

|/x|=3 



3! 



atl-.a^djX^iX) 



3! 



E -71— 7t(«'')'''^-(«''')'''''^^?>(^^) 



l/^'l=3 



- y 

2!1! ^ 



EE 

(=1 |/:t|=2 



2! 



2! 



+ ^ E 



1 



E E 

l',q' = l \tJ.'\ = l 



Mi!- 



where ^{^^ ( J2{i^j^k}) denotes the sum over ah different pairs (triples) of i,j G {1.2, .... f/} ( of i, j, k € 
{1, 2, d}) and G M'' denotes the vector whose i— th coordinate is equal to 1 and other coordinates are 
zero. Collecting the similar terms in the last equation we obtain that for v = 3efc, u' = 3ei the coefficient 

before x„{X)D';^' (f){AX) is equal to ^(aifca'^ + ... + adfca''')^ = j^Ski, for ly = eq + 2er,v' = e/ + 2e„ , g ^ r, 

the coefficient before Xv{^)^x 't'i-^^) is equal to ^{axqa''^ + ...+adqal''^){a\rCi'^^ +---+adrO'^'^)'^ = -^^^qi^m, 
in particular , for I = n the last expression is equal to zero. For = + + e„, i/' = e^' + e^' + e„' 
I'lQ ^ ''^I'l' ^ the coefficient before x^(X)_D^'(/)(j4X) is equal to (aiga^'-^ + ... + ad^a'''^) x {airOJ''^ + 
... + adra''''^)x (ai„a"'i + ... + adnO^'''') = dqq'drr'Snn' ■ This proves lemma for [ul — 3. The proof for 
\u\ = 4 is quite similar. For this case we use the relation which enabes to express a cumulant Xi^iAX) as 
^j^(ylA') plus a second order polynomial of the moments ^^,,{AX), \u'\ =2. A necessary correction term 
for Ht,{X) to get a Xvi^) comes from the derivation of D^(j){z). This completes the proof of the lemma. 

Lemma 5. The following bound holds with a constant C for v = (y\, ...Vp)^ with < < 6 

D'^Phi.jh, kh, x, y) - D'^pijh, kh, x, y) - VhD'^TTi{jh, kh, x, y) - hD'^n2{jh, kh, x, y) 
<Ch'/^p-\l-\'^\y-x) 
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for all j < k.xand y. Here D'^denotes the partial differential operator of order v with respect to 
z = Vj j}''^ {y){y — X — iiji^{y)). The quantity p denotes again the term p= [h{k — j)]^^'^ and the functions 
TTi and TT2 are defined in (4) and (5). We write Cp(') = P~'^C'^('/p) where 



C (^) = r,. . fci 



J[l + \\zt]-'dz' 



Proof of Lemma 5. We note first that Ph{jh, kh,x, •) is the density of the vector 



where, as above in the definition of the "frozen" Markov chain Yn, is a sequence of independent 

variables with densities q{ih,y,-), Pj kiv) — 12^=1 hm{ih,y). Let fh{') be the density of the normalized 



sum 



fe-i 



Clearly, we have 



Ph{jh, kh, X, •) = det [Vj^kiv)] fn{[Vj,k{y)] [■ - x - Pj^kiv)]}- 



We now argue that an Edgeworth expansion holds for fh- This impHes the following expansion for 
Ph{jh,kh,x,-) 

Ph{jh,kh,x,-) (9) 

S—S 

= det [VjMy]]-'^' [J^C" - 3r''Pr{-cl> : {x0,r})m,k{y)r"' [• - ^ - Hkiy)]} 



+[k-j\ 



r=0 



_ •i-(5-2)/2 



with standard notations, see Bhattacharya and Rao (1976), p. 53. In particular, Pr denotes a product of 
a standard normal density with a polynomial that has coefficients depending only on cumulants of order 
< r + 2. Expansion (9) follows from Theorem 19.3 in Bhattacharya and Rao (1976). This can be seen as 
in the proof of Lemma 3.7 in Konakov and Mammen (2000). 
It follows from (9) and Condition (A3) that 



\ph{jh,kh,x,y) - p{jh,kh,x,y) -h'^^\i{jh,kh,x,y) - hTT2{jh,kh,x,y) 

<Ch'l^p-'Cl-\^\y-x), 



(10) 
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where 



p{jh, kh,x,y) = det [Vj,k{y)r'^^ {2n)-P^^ 

exp{- J(y -X- ^ijJy)f [Vj,k{y)r^ [y-x- m,- fe(y))}, 



kl=3 



Tii{jh,kh,x,y) 
TT2{jh, kh,x, y) = p"^ det K-,fc(2/)]~^^^ 



^1 1^ X^^4 



where x^j,fe(2/) = s^j Xl.j.mC?')' = '^"^'^ ^""lul^iit °f PW3,k{y)] ^''^'i 



i+l,h 



{v—th cumulant of [y^-,fc(y)] ■Ci+i h}i ^-^d -^z'Piz) denotes the t^— th derivative of (j) with respect to 



— 1/2 

z = \Vj^k{y)] {y — x — Hj^kiy)) ■ it follows from the (conditional) independence of ~ j, k—1, 

that x.,j,k{y) = X xMX), where A = h'^^VjM'^^^ = = Cov{X,X),X = 



J2i=j ^i+i,h- By Lemma 4 for s = 3, 4 



E 



Dim = E ^^^;:v.A('^^/^-) 

\y\=s 

= (-l)V^ E ^^D'^MA{y -X- M,.,fc(y))) 



\y\=s 

= (-i)vE 



^.^'/'(K-.feCy)]"'^' (y - ^ - HkivM (11) 



where we put (^^(2) = ^(/i ^/^z), x,y(-^) = ^zj Ei=j Xvi^^^v)- It follows from (11) and the condition 
Bl that up to the error term in the right hand side of (10) the functions tti and 772 coincide with 
the functions tti and ^2 given at the beginning of Section 4. For ;/ = the statement of the lemma 
immediately follows from (10). For v > Q one proceeds similarly. See the remark at the end of the proof 
of Lemma 3.7 in Konakov and Mammen (2000). 

Lemma 6. Let L{d) be the set of symmetric matrices, and for < X~ < < 00 let 0^+ ^- C L{d) 
be the open subset of L{d) that contains all A € L{d) with \~ I < A < A"*"/. For A e L{d) define 
A = A{h) as the symmetric solution of the equation A? = A. Then for any k, I, i,j<d and A € ^- 
we have that with a constant Cm depending on m 

9" ay (A) 



{dXkiY 

Here aij (A) are the elements of A = A{A) . 



-(2m-l)/2 



(12) 



Proof of Lemma 6. For m = 1 the lemma was proved in Konakov and Molchanov (1984) (see Lemma 
4). Suppose now that (12) holds for m < I. From the equality AA = A we obtain for m = / + 1 



d^+\AA) = {d^+^A)A 



l + l 
1 



{d^A)dA + ... 



' + 1 
I 



dA{d^A) + A{d^+^A) = 0, 
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where d denotes elementwise differentiation of a matrix with respect to a fixed element of A. This implies 



{d^+^A)A + A{d^+^A) 



l + l 
1 



{d'A)dA 



l + l 
I J 



] dA{d^A). 



(13) 



Denote the symmetric matrix in the right hand side of (13) by A. Then equality (13) determines a linear 
operator £ mapping d'''^^A to A. In the linear space of symmetric d x d matrices we introduce the scalar 
product {X, Y) = trace(Ary). The operator £ determines a quadratic form 



{£X,X} = tTace[{XA + AX)X] = 2trace[X^X] > 2 V A" trace [ATX] = 2VX~{X,X}, 



where in the inequality we have used that A — v X I positive definite implies that X{A— v A I)X = 
XAX - V)rXX is positive definite. Similarly, we get {£X,X) < 2\fx^{X,X). Hence, 



\\£X\ 



< 2v/a+ 



and 



1 



< 



We obtain 



2V0^-" 
U+'A\\ < 

Using the induction hypothesis we get from (13) 



< 



2VA 



2VA- 



|A||. 



||d'+^^||<Q+i(A-) 



(2J+l)/2 



This completes the proof. 



From Lemmas 5 and 6 we get the following corollary. The statement of the next lemma is an extension 
of Lemma 3.7 in Mammen and Konakov (2000) where the result has been shown for < |6| < 2,a = 0. 

Lemma 7. The following bound holds: 

\D^DlpH{jh,kh,x,y)\ < Cp-I«l-I*lc^-W(y - a;) 

for all j < k, for all x and y and for all a,b with < |a| + |6| < 6. Here, p = [(fc — j)/i]^/^. The constant 
S has been defined in Assumption (A3). 

Proof of Lemma 7. For two matrices A and B with elements aij or h^i, respectively where aij{B) are 
smooth functions of bki we write \§§\ < C if 



dai 



< C for all 1 < i,j < d,l < k,l < d. To obtain the 

, 1/2 

assertion of the lemma we have to estimate the derivatives DyD^z, where z = f, {y){y — x — fj.j f^{y)). 

]^/2 1/2 1/2 

Note that z = z{V-^y. , iJ.j f,, x, y), where Vj ,. = V- ^, (y) and = ^j^kiv)- / = 1, 6 it follows 
from condition (Bl) and (8) that 



5V,.fc(2/) 



(dyy 



<Cp\ 



It follows from Lemma 6 that 



{dy) 

< Cp-(2'-l)/2 



<Cp^ 



(14) 



(15) 
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From inequalities (3.16) in Konakov and Mammen (2000) and from the representation of an inverse matrix 
in terms of cofactors divided by the determinant we obtain that 



< Cp 



-0+1) 



From (14)-(16) and from the chain rule we get 

""Ik 



{dyy 



Now, Lemma 5 implies the assertion of Lemma 7. 



< Cp-K 



(16) 



(17) 



4.2 Bounds on operator kernels used in the parametrix expansions. 

In this Section we will present bounds for operator kernels appearing in the expansions based on the 
parametrix method. In Lemma 8 we compare the infinitesimal operators Lh and with the differential 
operators L and L. We give an approximation for the error if, in the definition of Hh = {Lh — Lh)ph, the 
terms Lh and Lh are replaced by L or L, respectively. We show that this term can be approximated by 
Kh + Mh, where Kh = {L — L)ph and where Mh is defined in Remark 5 after Lemma 8 . The bounds 
obtained in Lemma 9 will be used in the proof of our theorem to show that in the expansion of ph the 
terms Ph '^h Hj['' can be replaced by ph 'S'h {Kh + M/i) 

Lemma 8. The following bound holds with a constant C 

\Hh{jh,kh,x,y) - Kh{jh, kh,x,y) - Mh{jh,kh,x,y) - Rh{jh,kh,x,y)\ 

<Ch^'^p-'C%y-x) 

with (p as in Lemma 5 for all j < k, x and y. For j<k — 1 we define 

K'h{jh, kh, X, y) = {L- L)\{x), Mh{jh, kh, x, y) 

= Mh,i{jh,kh,x,y) + Mh,2{jh,kh,x,y) + Mh^ri{jh,kh,x,y), 

Mh,i{jh,kh,x,y) = h^'^ ^^^{xAjh,x) - xAjhv)), 
Mh,2{jh,kh,x,y) = h ^ \; {xAiKx) - xA3h,y)), 



kl=4 



M'h^,{jh,kh,x,y) = -{Ll- L^)X{x), 



Rh{jh,kh,x,y) = h^^^ 22 — YM'n^r{jh,x)p^_^^{jh,x) - mr{jh,y)p^_^^{jh,y)] 

\v\=i ' r=l 



+5 J2 ^II("ifeO'/i.a;)-mfc(j/i,y))(i/fc / q{jh,x,9)h''-'"' {9) 



|i/|=5 k=l 



{l-ufD''\{x + uh{9))d.u 
Dl\{x 



d9+ / q{jh,x,9)h''{9) 



(1 - ufuD"+'"'\{x + uh{9))dv 



d9 



kl=4 



Y v\N{v,v')[m''' {jh,x)p^_^,{jh,x) - m"' {jh,y)iJ.^_^,{jh,y)\. 

k'l=2 
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Here is defined as L but with the coefficients "frozen" at the point x, Cj- denotes a d- dimensional vector 
with the r-th element equal to 1 and with all other elements equal to 0. Furthermore, for = 4, = 2 
we define 

N{u,u') = 2x[-'!=il+x[(— ')!=il-2^ 

where x(') is the indicator function. We put m{xY = mi{xY^ ■ ... ■rn,,i{xY'' and m,{xY = 0, j^! = 0. We 
define ii^{t, x) = J z'^q{t, x, z)dz and ii.^{t, x) = if at least one of the coordinates of v = [vi, Vd) is 
negative. We use also the following definitions 

A(a;) = ph{{j + '^)h,kh,x,y), 

1/2 

Here again p denotes the term p = [h{k — j)] . For j = k — 1 we define 

Kh{jh,kh,x,y) = Rh{jh,kh,x,y) = Mh,2{jh,kh,x,y) = M'f^.^{jh,kh,x,y) = 

and 

Mh,i{jh, kh, X, y) = h'^'^+^^Z^ ^q {j/i, x, h'^l'^iy -x- m[jh, - q {j/i, y, h'^/^iy -x- m\jh, . 



Proof of Lemma 8. As in the proof of Lemma 3.9 in Konakov and Mammen (2000) we have 

Hh{jh,kh,x,y) = Hl{jh, kh,x,y) - Hl{jh,kh,x,y), 

where 

Hl{jh,kh,x,y) = h-' J q{jh,x,e)[x{x + h{e)) - x{x)]de, (is) 

Hl{jh,kh,x,y) = h-^ J q{jh,y,e)[X{x + h{e))-X{x)]de, (19) 

h{e) = m{jh,x)h + 6h}''^,h{e) = m{jh,y)h + 6h^''^. 

For [A(x + h{6)) — X{x)] and [X{x + h(0)) — X{x)] in (18), (19) we use now the Taylor expansion up to 
order 5 with remaining term in integral form. To pass from moments to cumulants we use the well known 
relations (see e.g. relation (6.11) on page 46 in Bhattacharya and Rao (1986)). After long but simple 
calculations we come to the conclusion of the lemma. 

Remark 5. We show now that the function K'i^{jh, kh, x, y) + M'f^^^{jh, kh, x, y) in Lemma 8 is equal to 
Kh{jh, kh, X, y) + ^{Ll- 2LL + L'^)X{x) + M'^,^{jh, kh, x, y) where 

■rrh'^{3h,y) , 



MlSh,kh,x,y) = -h^ Yl ^^y iL-L)D^^X{x) (20) 
-3 V / {l-5fd5 j q{jh,y,e)^^^{L-L)Di'X{x + 5h{e))d6. 



1/^1=3- 



Thus in Lemma 8 we can replace K'i^{jh, kh, x, y) + M'f^{jh, kh, x, y) by Kh{jh, kh, x, y) + Mh{jh, kh, x, y) 
where Kh{jh,kh,x,y) = {L -L)ph{jh,kh,x,y), Mh{jh,kh,x,y) = ^{Ll - 2LL + L^)X{x) + M'^,M'^ = 
Mh,i{jh, kh, X, y) + Mh^2{jh, kh, x, y) + M'^ .^{jh, kh, x, y) and 

max{|M^(j/i, kh, x,y)\, \Mh{jh, kh, x,y)\} < Cp^^C,p{y - x). 
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= kh — jh. To show this we note that 



Ph{jh,kh,x,y) = j q{jh,y,e)X{x + h{e))de, 
where h{9) = m{jh, y)h + h}/^9. From the Taylor expansion we get 

Phijh, kh, X, y) = X{x) + hL\{x) + h''Y^ '"''^■^|''^^ £)^A(a; 



|/i|=2 



+^12 f i^~^?dS f q{jh,y,( 



)^^^D^'\{x + dh{e))d0 



and, hence, 



From 



K'hijh, kh, X, y) = Kh{jh, kh, x, y) + {L - L)[X{x) - Ph{jh, kh, x, y)] 

= Kh{jh, kh, X, y) + h{L^ - LL)X{x) + M''{jh, kh, x, y). (21) 



hiL^ - LL)\{x) + M'^ -iijh, kh, x, y) = h{L^ - LL)X{x) + ^{Ll- L'^)X{x) 

= ^{Ll-2LL + L')X{x) 

and from the definitions of the operators L, L and and from the Lipschitz conditions on the coefficients 
m{t, x) and a{t, x) we obtain that 



'-{Ll - 211 + 1"^) X{x) 



< Chp-%{y - x). 



Analogously, we have 



\U.\=2 



'{jh,y) 



(L - L)Di'X{x) 



< Ch^p-^Cp{y - x), 



(22) 



(23) 



3 ^ / {l-5fd5[q{jh,y,e)^^!^^{L-L)D''X{x + 5h{e))de 



|m|=3 

Now (21)-(24) imply the assertion of this remark. 
Lemma 9. The following bound holds: 



<Ch^'^p-\Ay-x). (24) 



Y,Ph ®h {Kh + Mh + Rhf''> (0, T, X, y)-J2ph ®h {Kn + Mhf> (0, T, x, y) 
< C{e)hn-'l^+'C%{y - x), 



r=0 



r=0 



(25) 



where limej^o C{e) = +oo. 

Proof of Lemma 9. For r = 1 we will show that for any £ > 

\Ph Ofc {Kh + Mh + Rh){0, kh, x,y)-ph(E>h{Kh + Mh){0,kh,x,y)\ 

= \Ph ®fc Rh{0, kh, x,y)\< Ch^'^-'{kh)-^/^+'B{^,e)C^p {y - x), p^ = kh. (26) 
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Clearly, to estimate ph 'Sih Rh{^',T,x,y) it is enough to estimate 



fc-2 

/i = /i3/2 ^ ^ / ~^(o^ z) - fUh, y))Dlph{{3 + l)h, kh, z, y)dz 

A — f\ J 



3=0 

for zy, |z/| = 4, and 

fe-2 



h'^Y.^j Ph{0,jh,x,z){f{jh,z)-f{jh,y)) f q{jh,z,e)h''-^'{e) 

X / {l-ufDlX{z + uh{e))dudedz 
Jo 



for z^, = 5, 1 < A: < d. Here f{t,x) is a function whose first and second derivatives with respect to x 
are continuous and bounded uniformly in t and x. After integration by parts we obtain 



= _/j3/2 Y^^l Dl'ph{0,jh, X, z){f{jh, z) - fijh, y))D:-^'ph{{j + l)h, kh, z, y)dz 

3=0 ■' 

k-2 

Y,h Dl^ph{0,jh, X, z)Dl''f{jh, z)Dl-^^-^^ph{{j + l)h, kh, z, y)dz 

3=0 

k-2 . 

^ /I / Pfc(0,i/i, X, z)Dl'+''f{jh, z)Dl-'^-<''ph{{j + l)h, kh, z, y)dz 



3=0 

for 1 < Z, s < d. Hence, 

k—2 

\h\ < Ch'/' J2 ^ Vjfe(fcLi/i) ^^^^ ~ - Ch'^'~'ikh)-'^'^^B{le)C'^{y - x). (27) 
In the same way after integration by parts we get with 1 < /, s < d. 

k—2 1 

^2 = -h''J2''l {l-u)Uu J de{m{jh,y)h'/^ + er--' J Dl^ph{0,jh,x,z){f{jh,z)- f{jh,y)) 

k-2 „! 



xq{jh,z,e)D''-^'Ph{{j + l)h,kh,z + uh{e),y)dz + h'^y^h [ (1 - / de{m{jh,y)h^''^ + B) 

X J Dl^ [p^O,i/i, a;, z)D''f{jh, z)q{jh, z, e)]D--^^-^'ph{{3 + l)h, kh, z + uh{e),y)dz 
ft'^ft^ {l-ufdu j d0{m{jh,y)h'/^ + 0r-''phiO,jh,x,z){f{jh,z)- f{jh,y)) 



3=0 

xDl^qijh, z, e)D--^^ph{{j + l)h, kh, z + uh{0), y)dz. 
It follows from (28) that 

\h\ < Ch^'^-'{kh)-^'^+'B{\,e)C^^{y - x). (29) 
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Claim (26) follows now from (27) and (29). For r > 2 we use the identity 

Ph <S)h {Kh + Mh + Rh)^'~^ (0, T, X, y) - Ph (3h {Kh + M^)M (0, T, x, y) 

= Ph <^h {Kh + Mh + Rh)^''-'^ - Ph <^h {Kh + M^)('--i)] (3h {Kh + Mh){0,T,x,y) 

+Ph^h {Kh + Mh + RhY^~'^^ ^hRh{0,T,x,y) 
= 1 + 11. (30) 



For r = 2 we obtain from (26) and simple estimate \{Kh + Mh){jh,kh,z,y)\ < Cp2^Cp^{y~ ^)^P2 
kh — jh, 



Oh 0/» {Kh + Mh + Rh) - Ph ®h {Kh + Mh)] ®h {Kh + Mh){Q, kh, x, y)\ 



fe-2 



<C'h'/^-'B{-,e)^h{jh)-'/'+%kh-jhr'^' / Ci{z-xKl{y-z)dz 



j=o 



< C^h'/^-^B{le)Bi^,e + i)(fc/^)^Cp (y - x) 



with p'^ = kh. For r > 3 we obtain by induction 



Ph Ofc {Kh + Mh + Rh) 



(r-l) 



Ph ^h {Kh + Mfc)(''-i) Oft {Kh + Mh){0, kh, x, y) 



< C^h^l^-^B{\,e)B{\,e + \)...B{\,e + !_I)(A/j)-+(-2)/2^S(^ _ ^) 



< r(e)/i3/2-^ 



[cr(i/2)] 



r(£ + 



e+(r-2)/2^S/ _ 



C;{y-x) 



with p = kh. To estimate // we use the following estimates 

\DlDlph{jh,kh,x,y)\<Cp-\-\-\'\C^^-\''\{y-x),\Dlph{^ 
\Dl{Kh + Mh + Rh){]h,kh,x + v,x)\ < Cp-\f{v). 



(31) 



(32) 
(33) 



The inequalities (32) and (33) are obtained by using the same arguments as is the proof of Lemma 7. 
Using these inequalities and mimicking the proof of Theorem 2.3 in Konakov and Mammen (2002) we 
obtain the following bounds for r = 0, 1, ... 

DlD'^ph ®h {Kh + Mh + RhY^'Ho, kh, X, y) 



< C^kh)-\^\-\''\+^B{^, i)B(l, i)...i?(^, i)C^I''l(2/ - X) 



< 



(34) 



Inoquality (34) allows us to estimate // = [ph 0/j {Kh + M'^ + i?/j)(''~^)] (g)/j Rh{0, kh, x, y). For this it is 
enough to estimate 



h^/^ Y^hj^h^h {Kh + Mh + Rh)^'-^^m3h, X, z) 
xD^Ph{{j + l)h, kh, z, y){f{jh, z) - f{jh, y))dz 



(35) 
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for r > 2 ,|!/| = 4, and 

n-2 

Y^h \ph(^h {Kh + Mh + Rhi'-'^^m Jh, X, z){f{jh, z) - f{jh, y)) 

X I q{jh,z,e)h''-^'{e) I {l-ufD''ph{{j + l)h,kh,z + uh{e),y)dudedz (36) 



for r > 2, |z/| = 5, 1 < ? < d. Here f{f,x) is a function whose first and second derivatives with respect 
to X are continuous and bounded uniformly in t and x. The upper bound for (35) follows from (34) by 
integration by parts exactly in the same way as it was done to obtain the upper bound for 7i, see (27). 
This gives the estimate 

h j[ph ®h {Kh + Mh + Rh)^'^-^'>]{0,jh,x, z) 
xD'Phiij + kh,z,y){f{3h,z) - fUh,y))dz 

< ne)h'''-' '^rff f {khr^^-^yx{y - x). (37) 

The upper bound for (36) follows from (34) by integration by parts in the same way as it was done to 
obtain an upper bound for I2, see (29). This gives for (36) the same estimate as in (37) and, hence, 

\II\ < cmh^/'-^^^^^^ikhr^^-'^Xiy - X). (38) 
The assertion of the lemma follows now from (26), (30), (31) and (38). 

Lemma 10. Let A{s,t,x,y), B{s,t,x,y),C{s,t,x,y)be some functions with absolute value less than 
C{t — s)~^/'^C,'\^{y — x) for a constant C and an integer k > S'd. Then 



^ A ® ^ (S + C) W {ih, jh,x,y)-J2^^h B^''^ {ih, jh, X, y) 

r=0 r=0 

00 

= Y,[A<S)h ^] <S)h [C<S)h *]^''^ iih,jh,x,y), 

where «' = E^o^^''^- 

Proof of Lemma 10. Under the conditions of the lemma all series are absolutely convergent. The 
assertion of this lemma is a consequence of the linearity of the operation 'E)h and of the possibility to 
permutate the terms in absolutely convergent series. 



5 Proof of Theorem 1. 

We now come to the proof of Theorem 1. Main tools for the proof have been given in Subsections 3.1, 
3.2, 4.1 and 4.2. From Lemmas 1 and 2 we get that 

n 

p{0, T,x,y) = Y,P® H^-^^ (0> T, x, y) + o{h''T)(j>c^^{y - x). 
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With Lemma 3 this gives 

p{Q, T, X, y) - Ph{Q, T, x,y) = + ... + Tr + o{h^T)^c,VTiy ' (39) 

where 



Ti = ^p®H('-)(0,T,ar,y)-5^p®ftifM(0,T,a;,y), 

n 71 

T2 = Y.P®hH^''\^,T,x,y)-Y,P®h{H + M'^ + VhN^)^-\Q,T,x,y), 

r=0 r=0 

n n 

n = Y.P®h{H + M'^ + VhN^f^\Q,T,x,y)-Y,P®h{H + Mh + y/hNii'''\Q,T,x,y), 

r=0 r=0 

n n 

Ti = Y,P'®h{H + Mh + VhN,)^-\0,T,x,y)-J2p'®h{Kh + MhY-\0,T,x,y), 

r=0 r=0 
n n 

n = J2P'^h{Kh + MhY''\0,T,x,y)-J2Ph<^h{Kh + Mh)^''\0,T,x,y), 

r=0 r=0 
n n 

Te = Y.Ph®h{Kh + Mht\Q,T,x,y)-Y,Ph®h{Kh + Mh + RhtH<i,T,x,y), 

r^O r=0 

n n 

TV = J2P'^®^^(^f^ + ^h + Rh)^''H0,T,x,y)-Y,Ph®hHl;;\0,T,x,y). 



Here we put Ni{s,t,x,y) = {L — L)7ri{s,t,x,y). 

We now discuss the asymptotic behaviour of the terms Ti, T7. 

Asymptotic treatment of the term T\. 

We start from the recurrence relations for r = 1, 2, 3, ... 

{<d,jh,x,y)-(p®hH^'^) {0,jh,x,y) 

= (p^iJ^'^-i)) ®hir] iO,jh,x,y) 

+ -(p(^hH^--^^)](^hH{0,jh,x,y). (40) 

By summing up the identities in (40) from r = 1 to 00 and by using the Unearity of the operations (S) and 
(g)/, we get 

{P - P"^) (0, jh, x,y) = {p(g)H -p(g)hH) (0, jh, x, y) 

+{p-p'')<S)hH{0,jh,x,y), (41) 

where we put 

DC' 

p'^iih, i'h, X, y) = ^(p ®h H^^'^ih, i'h, X, y). (42) 

r=0 

By iterative application of (41) we obtain 

{P - /) (0, jh, x,y) = {p®H -p(S)hH) (0, jh, x, y) 
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+ (p O - p Oh i?) 0/, $ (0, jh, X, y) , 

where i'h, z, z') = H{ih, i'h, z, z') + H H{ih, i'h, z, z') + ... = Y.'^^i H^''\ih, i'h, 
By application of a Taylor expansion we get 

{jp®H -p®hH) (0, jh, X, z) 

J-l i-(i+l)h 



= Y, [X{u)-X {ih)] dv 

j_Q Jih JRi 

= 2_. / {u— ih)du I X'{ih)dv 

+ V / {u-ihf [1-5] X" (s) \s=s, dvdSdu, 

,_n Jih Jo Jr'^ 



where A [u) = p(0, u, x, v)H{u,jh, v, z), Si = Si{u, i, S, h) = ih + 6{u — ih). 
Note that ^ 

/ X'{ih)dv= / —p{0,s,x,v)\s=ihH{ih,jh,v,z)dv 
J fid J fid as 

+ [ p{0,ih,x,v)-^H{s,jh,v,z)\s^ihdv= [ L*p{0,ih,x,v) 

jRd OS Jjid 

x{L — L)p{ih,jh,v, z)dv — / p{0,ih,x,v)[{L — L)Lp{ih,jh,v,z) 

Jr'^ 

—Hi{ih,jh,v,z)]dv= / p{0,ih,x,v)Hi{ih,jh,v,z)dv 

Jr"! 

+ / p{Q,ih,x,v){L'^ — 2LL + L'^)p{ih,jh,v,z)dv, 
where Hi{s,t,v, z) is defined below in (53). We get from (45) 



2_] I {u — ih)du / X'{ih)dv = —{pi^h Hi){0,jh,x,z) 
~^ Jih Jr'^ 2 



+ ^{j)®h Aa){Q,jh,x,z), 



where AQ{s,jh,v,z) = (i^ — 2LL + L'^)p{s,jh,v,z). The direct calculation shows that 

d 

4 



1 

Ao{s,jh,v,z) = - ^ {apg{s,v) - apg{s,z)){ari{s,v) - ari{s,z)) 



p,q,r,l=l 



"" IvM^dt ^ i (a,,(5,^;)-<7,,(.,.))(m,(5,t;)-m,(s,.)) 



d^p{s,jh,v,z) ^ 1 ^ dari{s,v) d^p{s,jh,v,z) ^ 

dvj)dv„dvr 2 ' dv„dvrdvi ~ ' 
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where we denote by (< 2) the sum of terms containing the derivatives of p{s, jh,v, z) of the order less or 
equal than 2. Note that for a constant C < oo and any < e < 5 



h 



(pO/i Hi){0,jh,x,z) 



< Ch(f>c,^ (z-x), 



-^{p<SihAo){0,jh,x,z) 



<C{e)h'/'j-^'/'-^Uc,VJKi^-^)- 



(48) 



First inequality (48) follows from (Bl) and the well know estimates for the diffusion density p and for 
the kernel Hi . The second inequality (48) follows from (Bl), (47) and the following estimate 



h f 

-^y^h / p{0,ih,x,v 



d^p{ih,jh,v,z) 



dv 



< 



< 



i=0 

2 



d^p{0,jh,x,z) 



dVqdVrdvi 

d^p{0,jh,x,z) 



dVqdVrdvi 



+ 



dp{0, ih, X, v) d p{ih, jh,v, z) 

dVn 



dVqdVrdvi 



dvrdvi 

+ Ch'/^j-a/^-^)B{^,e)ct>c,^iz-x). 



dv 



(49) 



Now we shall estimate the second summand in the right hand side of (44). Clearly 

A"(s) = ■^p{0,s,x,v)H{s,jh,v,z) + 2—p{0,s,x,v) 

>^-Q^H{s,jh,v,z)+p{0,s,x,v)-^H{s,jh,v,z). (50) 
Using forward and backward Kolmogorov equations we get from (50) after long but simple calculations 

J-l f{i+l)h 

u - ihf / (1 - <5) 



J2 / {u-ihf (l-S) A"(s)U=,, 



, dvdSdu 



i-i i.{i+i)h 



i=0 ''ih Jo fe^i JR'^ 



p{0,s,x,v)Ak{s,jh,v,z) |s=s. dvdSdu, 



(51) 



where 



and 



Ai{s,jh,v,z) = {L^ - 31"^! + SLL"^ - L^)p{s,jh,v, z), 
A2 = [LiH + 2LHi){s,jh,v,z), 

A3{s,jh,v,z) = [{L- L)Li +2{Li - Li)L]p{s,jh,v,z), 
A4{s,jh,v,z) = H2{s,jh,v,z). 

Hi{s, t, V, z) = {Li - Li)p{s, t, V, z) 

1 / d^'Oij (s, v) d^Oij (s, z) \ d^p{s, t, v, z) 

~ 2^ 



(52) 



dvidvj 



d'"mi{s,v) d''mi{s,z)\ dp{s,t,v,z) 



dvi 



1 = 1,2. 



(53) 
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Using integration by parts and the definition (52) of ^2, ^3 and A4 it is easy to get that for any < e < 1/2 
and for fc = 2, 3, 4 

JZ^ l'ii+l)h rl r- 

2_\ I (u — ih)^ / (1 — (5) / p{0,s,x,v)Ak{s,jh,v,z)\s=sidvd5du 

"^ih Jo JR'^ 



<Cie)h'/'-^<f>c,VJKi^-^)- 
For fc = 1 we shall prove the following estimate for any < £ < | 

\J / (u — ih)^ / (1 — 5) / p{0,s,x,v)Ai{s,jh,v,z)\s=sidvd6du 

Note that the function Ai{s, jh,v, z) can be written as the following sum 

1 

Ai{s,jh,v,z) = - ^ (^^^(s,?;) - o-ij(s,z))(crp5(s,i;) - (7j,5(s,z))((7zr(s,-y) 



(54) 



(55) 



-air{s,z)) 



i,j,p,q,l,r=l 

d^p(s,jh,V,z) 3 / / X / ^^/ / X 

+ 1 {aij{s,v)-aij{s,z)){apg{s,v) 



dvidvjdvpdvqdvidvr 4 . 



-a-pg(s, z)){mi{s, v) - mi{s, z)) 



d^p{s,jh,v,z) 



3 , dapg{s,v) 

2^ <^ij{s,v)- 



dvidvjdvndvndvi 4 ^ dv. 



(a,.(s,.) - c.,.(a,.)) /g.^'i^'::i + (< 4), 



dvjdvpdvqdvi dVr 



(56) 



where we denote by (< 4) the sum of terms containing the derivatives of p{s,jh,v,z) of the order less or 
equal than 4. By (Bl) and (56) it is cloar that the estimate for the left hand side of (54) for fc = 1 will 
be the same up to a constant as for the following sum for fixed p, q, r, I 



1=0 



ih 



[u-ihff\l-6)( p(0,.,.,.)^^^^4^ 

Jo jRd dVpOVqdvidVr 



8=8, dvdSdu 



After integration by parts w.r.t. Vp and with the substitution hw = (u — ih) in each integral we obtain 



3-1 l'{i+l)h 

{u-ihf I (1-6) I p{0,s,x,v) 

i-l i-ii+l)h 



d^p{s,jh,v,z) 



/ {u-ihf / (l-<5) / 

• n Jih Jo Jh 



fid dVpdVqdvidVr 

dp{0, s, X, v) d^p{s,jh, V, z) 



=s, dvdSdu 



Rd 



<C.Vc,^(.-x)/'»^/'(l-.)g.-^ 



dVqdvidVr 



1 



^ ^/ih + 6hw [{j - i)h - 5hw)fl'^ 



dvdSdu 



dddw 



'0 Jo 

r-l /.I 



Vih + 5hw [{j - Sw)h - ih)]^- 



-dSdw 
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<Chr^'/^-^^B{\,e)<l>c,V3E{z-x), 



X] I w'^dw [ {l-6)^/'^-^d5 I 
Jo Jo 



dt 



Vi[{j-l)h-t]^-^ 



(57) 



where B{p,q) is a Beta function and 4>c^p [z — x) is defined in Lemma 2. As we mentioned above (55) 
follows now from (57). By (B2), (44), (46), (47), (49), (54) and (55) we obtain for any < £ < i and 
j = 1,2,. ..n 

\{p<» H - p:»h H){0,T,x,z)\ < C{e)h^/^n-^^/^-'Uc,VTi^ - (58) 

We use now the following estimate for ^{ih,i'h,z,z') that was proved in Konakov and Mammen (2002) 
(formula (5.7) on page 284) 



mih,i'h,z,z')\<c^ 



i'h-ih ^) • 



From (B2), (44), (46), (57), (58) and (59) we obtain the following representation 

{p - p'^) {0,T,x,y) = ^{p^h Hi)iO,T,x,y) + ^{p^h Ao){0,T,x,y) 

+ ^(l"»h ffi ^) (0, T, X, y) + '^{p ®h Ao Oft $)(0, T, x, y) 
+R{0,T,x,y), 
where for any < £ < 1/2 

\R{0, T, x,y)\< C(£)(/i3/2- + /in-(i/2--))0^_^ _ x) 

= <t>c,VTiy-x)o{h'+'). 
This representation impHes that 

2^1 = ^\p0h{L'' -2LL + L^)p0h^)]{O,T,x,y) 

+ ^[p ®ft (L' - L')p ®h *](0, T, X, y) + Rt{0, T, x, y), 

where for any < £ < 1/2 

\RT{0,T,x,y)\ < C{e)hn-^/^+'cf>c^^{y-x) < C{e)h^+' (j>c^^{y - x) 



(59) 



(60) 



(61) 



(62) 



for (5 > small enough and where ^{s,t,x,y) = Y^^^QH'^^\s,t,x,y). Here the summand H^^\s,t,x,y) 
is introduced to shorten the notation. By definition we suppose that g i^h -ff {s, t, x, y) = g{s, t, x, y) for 
a function g. Note, that in the homogenous case (Jij{s,x) = aij{x),mi{s,x) = mi{x) and thus the second 
summand in (61) is equal to 0. 

Asymptotic treatment of the term T^. We will show that 

oo oo 

T2 - 3 0ft fW(0, T,x,y) + Y.P ®h + ^/^.i + T, x, y) 

r=0 r=0 

oo oo 

+ Y,P®h{H + Mh,2)^^\0, T, x,y) + Y,P<^h{H + M'^^s)^^\o, T, x, y) 

r=0 l-=0 
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with some positive S > 0. Note that it is enough to consider the case r > 2 because for r = 1,2 the 
left hand side of (63) is equal to zero. Note that (63) immediately follows from the following bounds for 



r = 2,3,, 



and 



p®h{H + M'^ + VhNi)^'^\0, T, X, y) 

-p (g)^ {H + Mh.i + Mh,2 + V)^^i)('"H0, T, X, y) 
-\p®^{H + Ml^Y^^ ~p®u i?^""^] (0, T, X, y) 

< Cie)h'^'-''^^T^^+^CvkHiv - X), 



p®h{H + Mh.i + Mh.2 + y^iVi)('-)(0, T, X, y) 
-p ®h (H + Mh,i + VhNi)^^'^ (0, T, X, y) 
-[p ®H {H + Mha)^'-'^ -p®h H^'-^i^, T, X, y) 



(64) 



(65) 



for all sufficiently small £ > with a constant C{e) that fulfills hm^ >o C{e) = +oo. First we prove 

the bound (64). Denote the expression under the sign of the absolute value in (64) by Tr- Note that 
To = Fi = 0. For r > 2 we make use of the following recurrence formula 

F^ = Tr-i(^hH+^p(^h{H + M'; + s/hNx)^''-^^ 

-p^h {H + Mh,i + Mh,2 + VhNxf''-^^] {M'l + \PhNx) 

+ [p ®h {H + Mh,i + Mh,2 + V^iVi)(^-i) -p(g>h{H + m;:'_3)('-i) 
= / + // + ///. 
We start with bounding // . First we will give an estimate for 

p®hiH + Ml + y)^Afi)('-i) -p®h{ii^ Mh,i + Mh,2 + VkNi)^"--^ 

For r = 2 we have p ^h^aiO, kh, x, y). It follows from (20) that it is enough to estimate 



0/, M'' 



k-2 „ 

Ji = h?^h / p{0,ih,x,v){f{ih,v) - f{ih,y))D!;;ph{{i + l)h,kh,v,y)dv 



(66) 
(67) 

(68) 



for |z/| = 4 and 



J2=/i'/'^/i j mih,x,v){f{ih,v)- f{ih,y)) j q{ih,v,e)9'' j\l-Sf 

x_Dj;+ei+e,^^((j + l)h, kh, V + Sh{e), y)d6dedv (69) 

for |!/| = 3. Here f{t,x) is a function with D^f{t,x),\v\ = 0,1,2,3 bounded uniformly in (t,x). An 
estimate for Ji follows from (27). This gives 



|Ji| < Ch'-%khr-'/'B{^,s)C%^{y-x). 



(70) 
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The estimate for J2 can be obtained analogously to the estimate of I2 (see (29)). By integrating by parts 
we get 



k '2 n\ p n 

h = h^/^Y,h J {1-Sfd5 J de-e"" J Dl'+''^\p{0,ih,x,v) 



The derivative 



X {f{ih, v) - f{ih, y))q{ih, V, 9)]D'^p^{{i + l)h, kh, v + Sh{0), y)dv. 



0^+"" [p{0. ih, X, v)ifiih, v) - f{ih, y))q{ih, v, I 



is a sum of 9 summands. By using integration by parts once more for all summands which contain 
D^p{0,ih,x,v) with < 2 , we obtain 



fe-2 



Yh ,.,l, X — — ^ < Ch^/^-^'B{e,e)ikhf'-\''^zl{y - x) 



(71) 



for any e G (0, 1/4). It follows from (70) and (71) that for r = 2 (67) does not exceed Ch^/'^-'^^B{e, e){khf-'^C,^^{y- 
x). For r > 3 we use the recurrence relation 

p^h{H + Ml + VhNi)^^-^^ -p(^h{H + Mh,i + Mh,2 + VhNi)^^-^^ 
= p^h{H + M'l + \/^Ari)(''-2) -p(^^^{H^ Mh,i + Mh,2 + VhNi)^^-^^ 



= I' + II'. 
From (72) we obtain for r = 3 



(72) 



k-2 



< Ch^/^-^'B{e,s)C^^{y-x)J2Hihy-Hkh-ih)-^/^ 



i=l 



< Ch?/^-^^B{e,e)B{\,e){khf-'l\'^{y-x). 
To estimate //' we use the following estimates 

D'iDliH + Mh,x + Mha + ^fhNi){jh,kh,x,v)\<Cp-^-\''\-\''\Cp{v - x), (73) 
Dl{H + Mh,i + Mh,2 + ^N^){jh,kh,x,x + v) <Cp-\p{v~x). (74) 

To prove (73) one can use the following estimates for the summands in Mh,i,Mh,2 and VhNi 

h^/^DlDl [Dlphiij + l)h,kh,x,v){f{jh,x) - f{jh,v))] \ < Cp-i-l«l-l*ICp(^;-x),|z/| =3, 

\hDlDl [D''^Ph{{j + l)h, kh, X, v){f{jh, x) - fijh, v))]\ < Cp-i-l«l-l''ICp(t' - x), 1^.1 = 4, 
h'^'D^Dl [D';+'^^+^^pu{{j + l)h.kh,x,v)p'{f{jh,x) - f{jh,v))] < Cp-i-l«l-l\,(z; - x), W\ = 3, 



for a function f{t. x) with \a\ + \h\ derivatives w.r.t. x that are uniformly bounded w.r.t. t. These estimates 
are direct consequences of Lemma 7. To prove (74) one can use the following estimates for the summands 
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in Mh,i,Mh,2 and VhNi 

h}'^Dl [Dlph{{3 + l)h,kh,x,y) |^=,+, {f{3h,x)-f{jh,x + v))]\<Cp-^Cp{v-x)M=2, 

\hDl[D^Jh{{j + l)h,kh,x,y) \y=,+, {f{jh,x) - f{jh,x + v))] \ < Cp-\p{v - x),\u\ = 4, 

h^l^Dl [Dl+^o+^-phiij + l)h, kh, X, v) p\f{jh, x) - f{jh, X + v))] I < Cp-\^{v - x), = 3. 

Again, these estimates follow from the estimates obtained in the proof of Lemma 7. Note that with 
4Vrk^'^{y)fij^k{y),x,y) = Vr^^^{y){y -x- fij^kiv)) it holds 



dz 
dy 



< 



C 



dz 
dx 



and that with z{y-^^^'^{x + v), IJ-j^k{x + v), x,x + v) = V-^^''^{x + v){v — p-j^ki^ + v)) it holds 



< 



1/2, 



c 



dz 



dx 



< 



dz 



dV, 



-1/2 



dV, 



-1/2 



dx 



+ 



dz 



dp 



dx 



<c(|HI + i). 



Now with the inequalities (73), (74) we can proceed like in the proof of Theorem 2.3 in Konakov and 
Mammen (2002). This gives the following estimate for r = 3,4, ... 



DlDl\p<^h {H + Mh,i + Mh,2 + VhNi)^^-^^]{jh, kh, X, v] 
< C^B{1, 1) X ... X bC-^, ^)/-'-l»l-l''ICp(t; - x). 



(75) 



Now we denote pi^r = P'^h {H + M^^i + Mh^2 + VhNi)'^'^\pQ = p. To estimate pi,r-2 ®h ^h^s it is enough 
to make the same calculations with integration by parts as it was done above for Ji and J2. This gives 



\II'\ < \pi,r-2<^hM'^^3{0,kh,X,y)\ 

< C"-/i3/2-2^B(l, i)i3(l ^ ^ 
and by induction 



(76) 



2' 2 



) X ... X B{1 



|/'| <C"-/i3/2-2-B(£,i)B(£+ ^ 



3 1. 



'■,e){kh) 



2' 2^ .... X B{s + -)B{e,e){khr+^CvkK{v - x), 

r = 3,4, .... Comparing (76) and (77) we obtain that for r > 3 

p^h{H + M'^ + \/^Ari)(''-i) - p 0^ (F + Mh,i + Mh,2 + VkNi)^"--^^ 



(77) 



< C^h^/^-'^B{e, l)B{e + ^, ^) x ... xB{e+'-^, ^)B{e,e){khr 



From (78) we get the following estimate for // 

|//| < C"-/i3/2-2-i?(£, e)B{e, l)B{e + ^, ^) x ... xB{e+'-^, ^)T'+^CM^ - x). 



r-2 1, 
2^ 



(79) 



2' ' 2' 2' ' 2 

To estimate /// note that the following inequalities that are similar to (73), (74), (75) hold for iJ + M^'g 

\D-Dl{H + M'^^^){jh,kh,x,v)\ < Cp-'-\-\-\^\Cp{v-x), 

\Dl{H + Ml^){jh, kh,x,x + v)\ < Cp-'Cp{v - x), 

D^Dllp^h {H + M'l^)^-)]{jh, kh, X, v) 



< C"-B(l, \) X ... X B(^, i)p'-H-l''lc,(^ - :.). 



(80) 
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To prove the last three inequahties it is enough to get the corresponding estimates for summands in M^' 3 
(see (20)). These estimates can be proved by the same arguments as used in the proofs of (73), (74), and 
(75). To estimate /// we have now to estimate pi^r -^^'3 and P2,r ®h ^h^s where 

P2,r=P®h (iJ + M;:'^3)W. 

Using integration by parts and inequahty (80), we obtain for P2,r-i ®h 3 the same estimate as for 

Pl,r (S^h K',3 



\p2,r M'^^{0, kh,x,y)\ 
< i) X ... X bC-^, i)B(^ 



for i = 1,2. Hence for r = 2, 3, ... 

\III\ < \pi,r-i'S)hM;i.^{0,kh,x,y)\ + \p2,r-i<»hM;i.^{0,kh,x,y)\ 

< C^h'/^-'^B{l, i) X ... X B(^, i)BC-^,e)ikhr+'^Cvkkiy - x). 

From (66), (79) and (81) we get for r = 2,3, ... 

|r.(0, kh, x,y)\< C^h'/^-''B{e, e)B{e, i) x .... xB{e + ^, lmy+'^Cvkk{v - x). 



(81) 



In particular. 



|r,(o, r, X, y) I < h^'^-'^ ^^^^ ^-+^Cv^(^^ - x), r = 2, 3... 



(82) 



for any e E (0, 1/4). Now we estimate the left hand side of (65). Denote the expression under the sign 
of the absolute value in (65) by Fr- Note that Fo = Fi = 0. For r > 2 we make use of the following 
recurrence formula 



r-l 



H 



p®h [H + Mh,i + Mh,2 + VhNif-^^ 



+ [p 0h {H + Mh,i + VhNi)^^-^^ - p<S)h {H + Mft,i)(^-i)] ®h Mh,2 
= I + 11 + III. 

We start again from the estimation of 

Ar-i = p<S)h {H + Mh,i + Mh,2 + VkNi)^""-^^ - p<S)h {H + Mh,i + VhNi)^''-^K 
For r = 2 we have Ai = {j)®h Mh,2){Q, kh, x, y). It is enough to estimate 
k-2 . 

J3 = h^h / p(0, ih, X, v){f{ih, v) - f{ih, y))D'^ph{{i + l)h, kh, v, y)dv 

for |!/| =4. Analogously to (27) we obtain that 

\Js\ < Ch'-'{kh)-y'+'Bi^,e)C^iy x). 
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and, hence, 



l^il < Ch'-%kh)-'/^+'B{^,e)C%j^{y - x). (83) 
For r > 3 we use the recurrence relation 

Ar-l = Ar-2 Oh {H + Mh,l + Mh,2 + VhNi) 

= I' + II'. (84) 
From (83) and (84) we obtain for r = 3 



fc-2 

|/'| < Ch^-'B{-,e)C%^{y-x)Y,h{ihf-^/\kh-ih)-^'^ 



i=0 



< Ch}-^B{\,e)B{\,e+\){khfe^{y-x). (85) 



^2' ' '2' 2 
To estimate //' we use the following inequality for r = 3,4, ... 

D'iDl^®h {H + Mh,i + VhN,)^--^^]{jh,kh,x,v) 
< C^B{1, 1) X ... X bC-^, l)p^-'-\'^\-\'\C,{v - X). (86) 

This inequality follows from (75). We have 

I J/'l < Ch}-'B{1, e){khyc'^{y - x). (87) 

Comparing (85) and (87) we obtain that 1^2] <C'^h}-^ B{\,e)B{\,e + \){khYC,^^{y - x). By induction 
we easily get for r = 2, 3... 



\Ar-^{Q,kh,x,y)\ < C-^h^-'BC-,e) x ... x B{^,e +'^)(khf+'^ (^^^(y - x). (88) 
Ar-i <S)h {Mh,i + Mh,2 + VhNi) 



it is enough to estimate 

fe-2 



J4 = /ii/2 f Ar-i{0, ih, X, v){f{ih, v) - fiih, y))Dlph{{i + l)K kh, v, y)dv 

i=0 



for \i>\ = 3, 

fe-2 



i=0 

for \v\ =4, and 

fe-2 



J5 = h^h / Ar-i{0,ih,x,v){f{ih,v) - f{ih,y))D'^ph{{i + l)h,kh,v,y)dv 



Jq = h^/^^h j Ar-i{0,ih,x,v){f{ih,v) - f{ih,y)){kh-ih)D^+''''+''''p{ih,kh,v,y)dv 

i=0 
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for If] =3. It follows from (88) that 



r-2 



r - 1 



£ + ^)B{e,e+ ■-—){khf^+'-^C.^{y - x). 



kh^ 



Clearly, the same estimate holds for J5 and Je . Thus we obtain 

\II\ < C^h'/'-'^B{l,e) X ... X B{l,e+'-^)Bie,s+'-^){khf^+'^C'vkh{y-x). 
Now we give an estimate for ///. We write 

Br-i = P®h {H + Mh,i + VhN^i''-^^ - p®h {H + Mh,it-^\ 
Using the recurrence equation 

Br-i = Br-2 (E>h {H + Mh,i + VhNi) +p(g,h{H + Mh,i)^''-'^^ (E>h VhNi, Bq = 
we obtain that 

r-2 

III = ^P3,i 0h VhN, {H + Mh,i + ^7Vi)('^-'-2) ®^ Mh,2{0,T,x,y), 

1=0 



(89) 



where psj = p®h {H + Mft,^i)('^. To estimate III it is enough to estimate a typical term in the last sum. 
Thus we have to estimate 

n—2 „ k—1 r r 

h"^Y.^ / Y.^ PsMih,x,w){jh-ih)D^+^-+^"^p{ih,jh,w,z) 

fc=0 '' j=0 '' i=0 '' 

X {g{ih, w) - g{ih, z))dw\ [H + Mh,i + VhNi)^'-^-'^\ih, kh, z, u)dz}(/(A:/i, v) - f{kh, y)) 
xD'^Ph{{k+l)h,T,v,y)dv. 

To estimate this term we apply two times an integration by parts in the internal integral / ...dw and then 
we make two times an integration by parts in / ...dv. We also use the following estimates 



DlD%,{0,ih,x,w)\ 



,1 + 1 1, ,., . 



< C'B{1, -) X ... X B{-^, ^){ih)^^C^{w - x), 
D^Dl{H + Mh,i + VhNr)^'-^-'^\ih, kh, z, v) 



,1 1 



< C"^-'-^B(- -) X ... X B{-, 



1 r-/-3 



){kh-ih) 



^-i-4-|„|-|6| 



^2'2' '2' 2 

for < Z < r — 3 with B{^, 0) = 1. This gives the following estimate for any 0<Z<r — 3, r>2 
m ®h ®h {H + Mh,i + \/hNi)^'-^-'^'^ 0h Mh,2{0,T,x,y) 

< C"'/,3/2-3e r(g) 3,+ r_3 _ 

r(3£ + i^) ^^^^'^ ' 



(90) 



For 1 = 1 — 2 we have to estimate 



P3,r-2 (®h V^iVlft Oft Mft,2(0, T, x, y). 
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This is a finite sum of terms corresponding to the different summands in Ni^h and Mh^2- To estimate a 
typical term 



ri-2 



k-1 



ft3/2^/^ / {E'^ fp3,r-2{0,jh,X,w){kh-jh)D!i;^^-+''-p{jh,kh,W,v) 
fc=0 j=0 

X {g{jh, w) - g{jh, v))dwYfikh, v) - f{kh, y))D-^Ph{{k + T, v, y)dv 



again we apply integration by parts and after direct calculations we obtain the following estimate for 



r = 2,3,... 



P3,r-2 O/i y/hNih 0ft M/,,2(0, T, X, y) 

The inequalities (64) and (65) follow now from (82), (90) and (91). 
Asymptotic treatment of the term T3. We will show that 



T3 



J2p^hiH + A)('-) {0,T,x,y)-J2p^h H^''^ (0, T, x, y) 



.r=0 



r=0 



< 



Chn \^{y-x), 



where A = M^' - = -f (LJ _ 2LL + L^)X{x). Write 

Cr = p®h{H + M'^ + VhN^t\Q,T,x,y) 
-p<S)h{H + Mh + \/^Ari) W (0, T, X, y) 
-\p^h{H + A)('-) - p ®ft if M] (0, T, a;, y). 

Similarly as in (66) we have the following recurrence relation 



Cr-l <»h H - 



p^iSih 



-p (^h {H + Mh + VhNi)^'^-^A <^h {M'^ + ^fhN^ 



+ 



= I + 11 + III. 



With the notation 



we get 



Dr-i =p<S)h{H + Mh + VkNi)^"-^^ -p<S)h{H + A)^^-^ 



Dr-i = Dr-2 (E)h {H + Mh + VhNi)+ph (^h {H + A)^'-^) (Mh-A + VhNi). 
Iterative application gives 

/// = Dr-l (^h A 
r-2 

= P4,; 0ft {Mh-A + VhNi) ®^ {H + Mh + VhNi)^^-'-^^ ®fc A{0, T, x, y), 



(91) 



(92) 



(93) 
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where P4_j = pi^h (H + A)^''h This sum can be estimated in exactly the same way as the sum in (89). 
This gives for r = 2, 3, ... 



\III\ < c(.)/.^/^-^^^;^r3^+^Cv^(t; 



X . 



(94) 



To estimate // we write 

Er-i =p^h{H + M'^ + V^iVi)('-i) -p^hiH + Mh + ^iVi)('-i). 
For r = 2 we have Ei =p^fiA and analogously to (70) 

\E,\ < Ch'-'{khr-'/'B{^,e)C'^{y - x). 
For r > 3 similarly as in (72) we use the recurrence relation 

E^_i = Er-2^h{H + M'f[ + VhNi) + \p^hiH + Mh + VhNi)^^-^^^hA 
= I' + 11'. 

The terms /' and //' have a similar structure as the corresponding terms in (84) and they can be estimated 
similarly. This gives the following estimates for r = 2, 3, ... 



,1 



r-2. 



\II\ = 



Er-i ®h {M'^ + V^iVi)(0, T, X, y) 



< C{e)h'/--^^,^^T'^+'^C^{v-x). 

The claim (92) follows from (93), (94) and the last two inequalities. 
Asymptotic treatment of the term T4. We will show that 



(95) 



T4 = Y.P®hH^''\^,T,x,y) 

r=l 

00 

-Y,P®h[H + hN2p (0, T, X, y) + Rl{x, y), 

r=l 

with N2{s,t, X, y) = [L — L)'jT2{s,t, x, y), |-R^(a;, y)\ < Chn~^(^{y — a;) for 5 > small enough and with 
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a constant C depending on 5. For the proof of (95) it suffices to show that for 5 small enough 

n 

Y,P®h{H + Mh + VhNi + /iiV2)(") (0, T, X, y) 



< 



-Y,P®h iKh + Mh)^-\0,T,x,y) 

r=l 

(jk 



< 



r=l 



Y,P®h{H + Mh + VhNi)^^^ (0, T, X, y) 

n 

■ J2p(^f^{H + Mh + VhNi + hN2)^'-^ (0, T, x, y) 

r=l 

~ n n 

c/in-^CvT(y-^)- 

Denote D^fi = and 

£'3,m(0,i/i,a;,2/) = ^p®,,(ii:,, + Mfc)W(0,i/i,a;,y) 

r=l 
m 

- 0^^^ + + + hN2)^^\0,jh, x, y). 

r=l 

Then (96) can be rewritten as 

|£>3,n(0,T,ar,t/)| < C/in-^C^(2/ - a;). 

We now make iterative use of 

-Ds.m = £'3,m-l 0^ {H + Mh + VhNi + hNi) + Qm-U 

for m = 1, 2, where 



gm{0,jh,x,y) = - 



{H-Kh + VhNi + hN2) (0, jh, X, y) 



.r=0 

Sh,m ®h {L - L)dh{0,jh,x,y) 



with 



goiO,jh,x,y) = -p®h (H - Kh + VhNi + hN2){0,jh,x,y), 

dh = Ph-P - VhTTl - hTT2, 
m 

Sh,m{0,ih,x,y) = ^p(g)h{Kh + Mh)^''\0,ih,x,y). 

r=0 
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Iterative application of (98) gives 

ri-l 

I)3,n(0, T, x,y) = Y^ 9r ®h {H + Mh + VhN^ + ftTVs)^"-'^-^) (0, T, x, y). 

r=0 

To prove (96) we will show that 

gr ®h {H + Mh + VhNi + /i7V2)("-'^-iH0, T, x, y) | (99) 
< f^^hn-'(:%{y-x). 

For this purpose we decompose the left handside of (99) into four terms 

a^i = J2 h f gr{0,ih,x,u){H + Mh + \/hNi + hN2)^''-''-'^^{ih,T,u,y)du, 

0<i<n/2 

ar-,2 = / / Sh,riO,kh,x,v){L — L)dh{kh,ih,v,u) 

n/2<i<n 0<fc<i/2 

x{H + Mh + VhNi + /lA^z)'""''"^^ {ih, T, u, y)dvdu, 
Or.s = ^ ^ J J {L^ — L^)Sh,r{0,kh,x,v)dh{kh,ih,v,u) 

n/2<i<n i/2<k<i-n^' 

x{H + Mh + VhNi + hN2f"'^"^^\ih,T,u,y)dvdu, 
ttr^A = j J {L^ — L^)Sii,r{0,kh,x,v)dii{kh,ih,v,u) 

n/2<i<n i-n^' <k<i-l 

x{H + Mh + VhNi + hN2)^''-''-^\ih,T,u,y)dvdu. 

Here and denote the adjoint operators of L and L, and 6' satisfies inequalities 2k < (5' < |(1 — xr), 
where x is defined in (B2). For the proof of (99) it suffices to show for Z = 1,2,3, 4 

Ki\ < hn-'C"-^B{l, 1) X ■■■ X i?( " ~ p \ l)C^(t/ - x) (100) 



< 



r(^) 



hn-'C%{y - x) 



for some 6 > 0. 

Proof of (100) for I = 2. Note that k < i/2,i > n/2 imply ih- kh> ^. The claim follows from 
the inequalities 

ma.x{\Kh{ih,jh,x,y)\,\Mh{ih,jh,x,y)\, VTiNi{ih,jh,x,y) (101) 

\hN2{ih,jh,x,y)\ , \H{ih,jh,x,y)\} 
< Cp~^( {y — x) with = jh — ih for < i < j < n, 
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\Sh,m{0,kh,x,v)\ < CC^^iv - x), 

{L - L)dh{kh, ih, v, u) < Ch^^^{ih - A;/i)-2(Jzl_(u - v) 
= 0(/.n-VW^)C5^(.-.), 
{H + Mh + VhNi + /lATs )("-'-- 1) {ih, T, u, y) 



(102) 
(103) 

(104) 



,1 1- 



,n-r-2 1. 



< 



r(2 



for n-r- 3 = -1,0,1, ...,n- 3 with =T-ih. We put S(i,0) = 1). Inequality (101) follows from the 
definitions of the functions Kh,.-.,H. Inequalities (102) and (104) can be proved by the same method 
as used in the proof of Theorem 2.3 in Konakov and Mammon (2002) (pp. 282 - 284). Inequality (103) 
follows from the inequality ih — kh> -j, Lemma 5 and the arguments used in the proof of Lemma 7. 
Proof of (100) for I = 3. Note that n/2 < i,k > i/2 imply kh > -j. We use the following inequalities 



\dh{kh, ih, v,u)\< Ch^/'^iih - fc/i)-3/2(Jzl_(u - v) 
{L^ - L^)ShAO, kh, X, v) < CT-\lzl{v - x). 



i/2<k<i-n^' 



i/2<k<i-n^ 
fih—n^ h 



{L^ — LF)Sh,r{^, kh, X, v)dh{kh, ih, v, u)dv 
1 



(105) 



(106) 



;\ih-khY/^^^^'^ 



< Ch^/'^T-^ 



du 



hji {ih - m)3/2 



C^{u -x)< Ch^/'T-'-^''r^^'-''y\^{u - X) 



<Chn-' C^{u-x), 



where S" = S' /2 - x > 0. Claim (100) for Z = 3 now follows from (106) and (104). 

Proof of (100) for / = 4. For ?: - n''' < A; < z - 1, n/2 < i we have ih > T/2, kh > T/3, {i - k) < n^' 
for sufficiently large n. The integral 



(L^ — L^)Sh,r{0, kh, X, v)ph{kh, ih, v, u)dv 



is a finite sum of integrals. We show how to estimate a typical term of this sum. The other terms can be 
estimated analogously. We consider for fixed j, I 



I 



d'^Sh,r{0,kh,x,v) 
dvjdvi 



{aji{kh,v) - aji{kh,u))h ''/^ 



(107) 



iq^'-^^ [kh, u, h-^''^{u -v-hJ2 milh, u))]dv 



l=k 



-I 



-\fhu 



d^ShA^^; kh, X, v) 
dvjdvi 

x[aji{kh,u* — Vhw) — aji{kh,u)]q^^~^\kh,u,w)dw. 
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where u* = u — hY^\_l.m{lh,u). Now using a Tailor expansion we obtain that the right hand side of 
(107) is equal to 



Note that 



d^Sh,r{0, kh,x,u*) 



kl=i 

+2hY: 

kl=2 



dvjdvi 



^ d'^Sh,r{0, kh, X, u* - S\/hw) 



dvjdvi 



Dyji{kh,u) 



[w + Vhj:izlm{lh,u)]'' 



»1 i-l 

/ D!^aji{kh,u — SVhw — 6h'y^m{lh,u))d6 
•^0 l^k 



''\kh,u,w)dw. 



i-l 



dvj dvi 



l=k 



-h 



d^Sh,riQ, kh, X, u* 



dvjdvi 



rripQh, u), 



l=k 



d'^Sh,r{0,kh,x,u*) 
dvjdvi 



{wp + \/h^^mp{lh,u)){wq + Vh^^mq{lh,u)) (108) 



l=k 



l=k 



X 1^ Dyji{kh, u)d5 + [D>ji(A;/i, u - SVhw 



—Sh m{lh, u)) — D'^aji {kh, u) 



l=k 



dS > g^' ''\kh,u,w)dw 



^ d^ShA^,kh,X,U*) r (i-fe)(fc;,^^^^)rf^ 

ovjdvi J 



^ ^2 d'^Sh,r{0,kh,x,u*) 
dvjdvi 



i-l 



^2 iTipilh, u) mq{lh, u) + R, 



l=k l=k 

where by (A3') we have for jo < {i — k) < ,w' = {i — k)~^/'^w 

d'^Sh,r{0,kh,x,u*) 



\R\ < Ch^/^ 



dvjdvi 



n^'i'^ \\w'\\ + o(ri/2^-i/2+5'-)y 



< ChC^iu - X){h^n^^' + l)T-3/2„-l/2+3<572 j ||y,/||3 



(109) 
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We obtain analogously 



i—l 

h / Wp{wq + Vh'y^mg{lh,u))D^''aji{kh,u) l 
J i=k Jo 



^ d^ShA'^, kh, X, u* - SVhw) 



xg^' '^\kh,u,w)dw 
daji{kh, u) d^Sh,r{0, kh, x, u*) 



WpWgq 



dvpdvjdvi 
'•'^''\kh,u,w)dw + R, 



dS 



where 



duq dvpdvjdvi 

\R\ < C/in-i/2(i-35'-3>.)^^^^ -x),l- 3(5' - 3>^ > 



and, for 1 - 36' - 2>i > 

^3/2 Z"^^, d^Sh,r{0^kh,x,u* -SVhw) ^^ 
J Jo dvpdvjdvi 

'■^ d'^aji{kh,u — SVhw — 6hY^\Zjf.mq{lh,u)) 



durdi 



dS 



x{'Wr + Vh''^mr{lh,u)){ws + Vh'^^ma{lh,u))q^^ ''\kh,u,w)dw 

l=k l=k 

For 1 < i — k < jo the same estimates remain true because the following bound holds 

J \\wfq^^\t,x,w)dw<C{jo). 



(110) 



(111) 



The same estimates hold for p{kh,ih,v,u) with ''\kh,u,w) instead of q^^ {kh,u,w), where 
(j){kh,u,w) is a gaussian density with the mean and with the covariance matrix equal to a{kh,u). 
The first two moments of and (f)^^~^^ coinside so after substraction we obtain uniformly for 

i-n^' < k < i - 1 



{L^ - L'^)Sh,r{0,kh,x,v) 



n/2<i<n i—n^'<k<i—l 

X {ph{kh, ih, V, u) — p{kh, ih, v, u))dv 

x{H + Mh + VhNi + /lATs )("-''- 1) {ih, T, u, y)du 



(112) 



< 



r(^) 



/^y3/2^-3/2(l-.-55'/3)^^(^_^)_ 
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To estimate the other terms in dh{kh,ih,v,u) we need bounds for the following expressions 

h ^ I'iL'^ -L^)Sh,ri^,kh,x,v)Vhiih-kh) 



i—n^ <k<i—l 

xD'^p{kh,ih,v,u)dv for \u\ = 3, 



h J2 f{L^ -L'^)Sh,r{^,kh,x,v)h{ih-kh) 



i-n^' <k<i-l ' 

xD'^p{kh,ih,v,u)dv for = 4, 



(L' - L')Sh,r{^,kh,x,v)h{ih-khy 



i-n' <k<i-l 



We have 



X D^p{kh, ih, V, u)dv for \u\ =6. 



/i ^ /(^^ - L^)ShA'^^ kh, X, v)Vh{ih - kh) 



(113) 



D!^p{kh, ih, V, u)dv\ 



^ ^ [ ~ i'^)5'h,r(0, yfc/i, a;, - kh) 



i—n^' <k<i—l 



i— <k<i—l 



y/ih — kh 



<Chn-^^-^-^^'/^\^{u-x). 

Clearly, the same estimate (113) holds for = 4 and = 6. Now (100) for i = 4 follows from this 
remark and (112) and (113). 

Proof of (100) for I = 1. Note that for this case T-ih> T/2. 



ar,i= J2 12 {L^ -L'^)ShAO,kh,x,v) 

0<i<n/2 0</s<i-l"^ 

xdh{kh, ih, V, u)^h,r{ih, T, u, y)dvdu 



(114) 



0<k<n/2-l 



h jiL"^ - L'^)Sh,r{^,kh,x,v) 

h J dh{kh,ih,v,u)'^hAih,T,u,y)du 

^k+l<i<k+n^' 

+ h J dh{kh,ih,v,u)'i'h,riih,T,u,y)du\ 

k+n^' <i<n/2 ) 



dv. 



where we denote 



^hA^h,T,u,y) = iH + Mh + VhN,+hN2f'-''-^\ih,T,u,y). 
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We consider 



/i—l 
u — V — h'^^m{lh, u)])'^ h,r{ih, T, u, y)du 

J2 hj\ q^'-^^ {kh, v,w) + VhJ2('^ + '^12 "^(^^' u)YDlq^^-^'^ {kh, v, w) 

k+l<i<k+n'^' [ kl = l l=k 

kl=2 

^3^3/2 {w + y/hY.\Zlm{lh,u)Y 
kl=3 

X / {l-5fD''^q^^-''\kh,v + 5h^'^w + 6hy2m{lh,u),w)d6\ 
Jo l=k J 

f 

y \v\=i i=k 

kl=2 ^' 
kl=3 

X / (l-5)2l?;;'J'ft_^(i/i,T,t; + 5/i^/2w + (5/iVTO(Z/i,u),y)d<5lrfu; 

This integral is a sum of 4 x 4 = 16 integrals. We estimate only two of them. Other integrals can be 
estimated by similar methods. First, we estimate 

h j q^^~''\kh,v,w)'^hA'i'h^T,v,y)dw = ^ h^h,r{ih,T,v,y)dw. 

k+l<i<k+n^' k+l<i<k+n^' 

Note that we get the same term when we replace (jr('~'^) [kh, v, w) by {kh, v, w). After the replacement 
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this term disappears. Second, we estimate 



/i— 1 

g^*-*^) {kh, V, w)Vh + Vh^m{lh, u)f D^^-^hAi^, T, v, y)d' 

d „ i—1 

J=l fc+l<i<fe+ra''' '=fe 



fc+l<i<fe+ 

+0{hn^' \\w\\+h^''^n'^^')]dw 

d i-l 
i=l fe+l<i</c+n'*' 

\ J=l fe+l<i<fe+n«' 



+0 



\ J=l fe+l<i<fe+n«' 



tz;) llwll dw 



J = l fe+l<i<fe+Ti''' 



l=k 



where 



The first term in the right hand side of this equation will be the same if we replace g*-' ''\kh,v,w) by 
(j)^^~^\kh,v,w). After the replacement this term disappears. For a proof of this equation we consider 
the function u{w) that is defined as an implicit function and we used the following change of variables 



h^^^w = u — V — h m{lh, u) 



l=k 



to obtain 



i-l 



i-l 



V/^^to(I/i,wH) = Vh^m{lh,v) + O (h{i - k) \\w\\ + h!^'^{i - kf^ 

l=k l=k 

because of {i — k) < . By similar methods we get 

h I [VhTri{kh,ih,v,u) + h'iT2{kh,ih,v,u)]^h,r{'ih,T,u,y)du 



k+l<i<k+n^ 



(115) 



It remains to estimate 



h I dh{kh,ih,v,u)^hA^^j'^->'^->y)'^'^- 



k+n^' <i<n/2 
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From (104) and (105) we obtain 







<n/2 


< 


^n—r 


r("-^-') 


< 


^n—r 




< 




r("-^-i) 



(116) 



7^-l/2/j3/2 



T/2 



T-y^hn-'''\,j^{y-v) 



Now we substitute the estimate (116) into (114). This gives the following estimate for any < e < x 

2/2-1 



h / dh{kh,ih,v,u)'^h,r{ih,T,u,y)du 

' ' '2 

n/2 
fc=l 



(117) 



k+n^' <i<n/2 
^n—r 



n 



n—r—l > 



<C{e) 



r(^). 



For fc = we get with 5/(,r-(0, 0, x, v) = 6{x — v) where 5{-) is the Dirac function that 

j j Sh,r{0,0,x,v){L-L)dh{0,ih,v,u)'^h,r{^h,T,u,y)di 
hn-('/^-'\^{y-x). 



l<i<i/2 

<C{e) 



^n—r 



r(^) 



This completes the proof (100) for 1=1. The estimate (97) may be proved by the same arguments as 
were used in the treatment of T^. 

Asymptotic treatment of the term T5. We will show that, 

00 

T5 = -v^^5?i«)^(fl' + M^,i + \/^Afi)M(0,r,a;,y) 

00 

-hY,^2®h H^''^ (0, T, X, y) + Rh{x, y) , 



(118) 



r=0 



where \Rh{x,y)\ < Chn ^C^^(y — x) for some 7 > 0. Note that with Sh(s,t,x,y) = X)r=i(-^^ 
Mh)^^\s,t,x,y) the term T5 can be rewritten as 



T5 = {P- Ph){0, T, X, y) + (p-ph) ®h Sh{0, T, x, y). 
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We start by showing that for n < 5 < i-j^ uniformly for x,y G R 



h X! (.Ph-p){0,jh,x,u)Sh{jh,T,u,y)du 



1<J<" 



< 0{hn- 



_l/2(l-^-45)^.S-2/ _ 



)C^\y-x) (119) 



for S small enough. For the proof of (119) we will show that uniformly for 1 < j < and for x,y £ R'^ 
j Ph{0,jh,x,u)Sh{jh,T,u,y)du = Sh{jh,T,x,y) 

j p{0,jh,x,u)Shijh,T,u,y)du = Shijh,T,x,y) 

+0[/ll/2r-V2n-V2+5 + /jl/2j.-l + ^V2/jl/2]^S^2(^ _ (^21) 

Claim (119) immediately follows from (120)-(121). For the proof we will make use of the fact that for all 
^ < j < and for all x,y G R''^ and \p\ = 1 



ID^Shijh, T, X, y)\ < C{T - jh)-\^^^{y - x) 



(122) 



Claim (122) can be shown with the same arguments as in the proof of (5.7) in Konakov and Mammen 
(2002). Note that the function $ in that paper has a similar structure as Sh- For 1 < j < the bound 
(122) immediately implies for a constant C 



\D:Sn{jh,T,x,y)\ < C'T-\^^\y-x). 



(123) 



We havep,,(0,i/i,a;,u) = /i"''/2gW[0,u,/i-V2( 

u — X — h^l^Q m{ih,u))]. Denote the determinant of the 
Jacobian matrix of u — hJ2iZo''^{'^^''^) by ^h- From the condition (A3) and (123) we get that for 

1 < j < n'^ 

j Ph{0,jh,x,u)Sh{jh,T,u,y)du 

= / h-'^/^q^^^[0,u,h-^^'^{u - x - hY,m{ih,u))]Sh{jh,T,u,y)du 

■' i=0 

= / q'^'\0,x + h^/'^w + h^m{ih,u{w)),w) \ A~^\ Sh{jh,T,x + h^^'^w + h^m{ih,u{w)),y)dw 

= J [q^^\0,x,w)+O{j-''^\'/^){\\w\\+l)i;ir'^'w)][l + O{jh)^^^^^ 
+0(/.V2T-i)Cjz2(y - x)il + \\wf~^){\\w\\ + l)]dw 

= S^ijh, T, X, y) + 0[/i1/2t- V2„-l/2+^ + /,l/2y-l ^ ;^l/2^V2]^S-2(y _ ^) 

with u = u{w) in X]i=o ™(*^'^) defined by the Inverse Function Theorem from the substitution w = 
h^^f'^{u — X — hYyi^o m{ih,u)). This proves (120). Claim (121) follows by similar arguments. From 
(119) we get that for 5 < ^ (with x defined as in (B2)) 

T5 = {p-Ph){0,T,x,y) + h J{p-ph){0,jh,x,u)Sh{jT,u,y)du + Rh{x,y) 

n^<j<n 
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with \Rh{x, y)\ < 0{hn ^/'^^^ ^ — a;) . We now make use of the expansion oiph — V given in 

Lemma 5. We have with p = {jKfl'^ > W^l'^n^l'^ 



n ,. 

h h^'/^p-^ j Cf{u-x)ShUh,T,u,y)di 



< Ch'^T-^'n-^' 



n „ 

' J2p~'^''' J \cU''-x)Sh{jh,T,u,y) 



du. 



where 5' < U{1 - 6)-^ ,2 6" = 5 + 256' - 25'. Now we get that 



(124) 



n „ 

h P'^^^'' / \cl{u- x)Sh{jh,T,u,y)\du <CB{5' ,l/2)T''-^'\^^{y - x) (125) 



for a constant C. This shows that for 5' > Q small enough 

Is = -[\//i7ri + ft,7r2](0,T, a;,y) 



-h Y [y^TTi + h'K2]{0,jh,x, u)Sh{jh,T,u,y)du + R'f^{x,y) 



n° <j<n 



with \R'j^{x,y)\ < 0{hn ' ^^'^^)C,^^{y — x) with a constant in O(-) depending on 5' . It follows from 
(119), (124) and (125) that 



Ts = - ^[VT^^fi + hn2\ ®h {Kh + Mft) W(0, T, x, y) + R'l,{x, y), 



(126) 



r=0 



where |i?^'(x,y)| < 0{hn-'^^"-'^/'^'>)('^^^{y - x). Now we apply Lemma 10 with A = VhiTi, B = H + 
Mh,i + VhNi, C = {Kh-H- VhNi) + (M^ - M^,i) to 



- Y ^'^1 '^h {Kh + MhY''^ (0, T, t/) + ^ Vhni {H + Mh,i + VhNi)^''^ (0, T, x, y) (127) 

T-=0 r=0 

and with ^ = /i^2, B = H,C = {Kh - H) + Mh to 

OO CO 

- ^ hn2 ®h {Kh + Mhf'> (0, T, x,y) + Y ®h H^'^Ho, T, x, y). (128) 

r=0 r=0 

The estimate (118) follows from (125), (127), (128), Lemma 10 and Lemma 5 . 
Asymptotic treatment of the term Tq. By application of Lemma 9 we get that 

\n\<C{s)hn-'/'+%'^{y~x). 



Asymptotic treatment of the term T7. From the recurrence relation for r = 2,3, ... 

Ph Oft {Kh + Mh + Rhi''^ (0, T, X, y) - ph (^h Hh^''^ (0, T, x, y) 

= Ph 87, {Kh + Mh + Rhi""-^^ - Ph ®h Hh^""-^^] ®h Hh{0, T, X, y) 

+ [pft ®h {Kh + Mh + i?ft)("-i) ®h {Kh + Mh + Rh- Hh)]{0, T,x,y) 
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and from Lemma 8 with r = 1 we get by similar arguments as in the proof of Lemma 9 that 

\T,\ < Ch^l^T-'/\^^{y -x)= Chn-'/'c'vriy ' 



Plugging in the asymptotic expansions of Ti, T7. We now plug the asymptotic expansions of Ti, T7 
into (39). Using Lemma 10, Theorem 2.1 in Konakov and Mammen (2002) we get 

Ph{0,T,x,y) -p{0,T,x,y) 

= Vh[ni+ / (g)h 3?i] (E>h ^{0, T, X, y) 

[i^2 + ^1 ®h $ ®h 3*1 + / ®h 3*2 + / ®h 3*3] ®h 4'(0, T, X, y) 
V ®h (3?i ®h'i^f^ (0,T,.T,y) 

+\p (^h (Ll - L^)p\Q, T, X, y) - ^-p ®h {L' - Z')/(0, T, x, y)] 
+0{h'+'C^{y-x)), (129) 

where 

00 

p'^{ih,i'h,x,y) = ^p®h H^''\ih,i'h,x,y), 

r=0 

5ii(s,t,a;,y) = Ni{s,t,x,y) + M-i_{s,t,x,y) - M-i_{s,t,x,y), 
^2{s,t,x,y) = N2{s,t,x,y) +Ui{s,t,x,y) - Ui{s,t,x,y), 

^3{s,t,x,y) = 2^ — ^ ^ '-Dlp{s,t,x,y), 

Mr{s,t,x,y) = ^k^^D^J{s,t,x,y), 

M,{s,t,x,y) = Yl ^^''^D^J{s,t,x,y), 
kl=3 ^' 

Ili{s,t,x,y) = ^ ^"^^^^^^ D'^wi{s,t,x,y), 

ni{s,t,x,y) = ^^''^D^,n,{s,t,x,y). 

Note that for the homogenous case and T = [0, 1] (129) coincides with formula (53) on page 623 in 
Konakov and Mammen (2005). 

Asymptotic replacement of p'^ by p. It follows from (42), (57) and (58) that 

\{p' - p){ih,jh,x,z)\ < C{e)h'-%jh - ihr-'/'(j>^y^-^{z - x) (130) 

for any < e < 1/2. Using (130) and making an integration by parts we can replace p"^ by p in (129). 
For example the operator — is an operator of order three. Applying integration by parts we get for 
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\v\ =3 

n-l 



Y,hj D:p{Q, ih, X, z)(p'' - p){ih, T, z, y)dz < C{s)h'-' ^ h-^ 

i=l i=l ^ ' 



By (B2) we have < x < 1 - 4£. This implies 



h 



p^h {Lt-L'){p^-p){Q,T,x,y) 



1 {ihf/^T-ihy/^-^ 

< C{e)h^'^-^'T^^-^'^B{e, e + \)<t>^{y - x). 



< C{e)h'+'4>^{y-x) 



for some < (5 < 1/2. The other terms in (129) containing p'' can be estimated analogously. Thus we get 
the following representation 

Ph{0,T,x,y) -p{0,T,x,y) 

= Vh[nj^+p(g)h 3?i] (g)h $(0, T, X, y) 

[1^2 +^i®h'^®h^i+p(S)h^2+P ®h 3*3] ®h *(0, T, X, y) 
+P®h (3?i ®h $)^'^ (0,r,a;,y) 

+\p ®h {Ll - L')p{0, T, X, y) - ®h {L' - L')p{0, T, x,y)] + 0{h'+'C^{y - x)). 

In the further analysis we make use of the following binary operation (g)^. This operator generalizes the 
binary operation (8) introduced in Konakov and Mammen (2005). For s € [0,t — h] and t € {h, 2h, T} 
the operation ig)^ is defined as follows 

s<jh<t—h 

Note that for s G {0, h, 2h, T} the two operations (g)^ and 0/i coincide. 

Asymptotic replacement of (p(8)ft3?i)(8)ft$(0, T, a;, y) 6»y p0(3?i0';j$)(O, T, x, y) = (p(8)3?i)(8)ft$(0, T, x, y), z = 
1, 2, 3, [p (3?i $)] (3fJi ®h $) (0, T, X, y)hyp® [(3?,J3'^ *) (3?, $)](0, T, a;,j/), p ®h {Ll - 
L^)p{0, T, X, y)hyp® {Ll - L^)p{0, T, x, y) and p ®h {L' - L')p{0, T, x, y) by p® {L' - L')p{0, T, x, y). 

These replacements follow from the definitions of 3?,, i = 1, 2, 3, and can be proved by the same method 
as in the treatment of Ti. There an estimate for the replacement error of p ®h H by p ® H is given. 
Linearity of the operation imphes that it is enough to consider the functions p®h^ where S>(w, t, z, v) 
is a function that has one of the following forms: 

^ktA^^^^J^D^,p{u,t,z,v) with |H=3,4, 

{L — L)tti{u, t, z, v) or {L — L)tt2{u^ t, z, v). 
We consider the case '^{u,t, z,v) = {L — L)TTi{u,t, z,v). The other cases can be treated similarly. It is 
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enough to consider a typical term of {L — L)'K\{u,t,z,v). We will give bounds for 

du [ p{0,u,x,z) ( [ xA'^,v)dw \ D'^{L-L)p{u,jh,z,v)dz 



J p{0,ih,x,z) y^J^^ xAw,v)dwj D'^{L - L)p{ihjh, z,v)dz 

i-i i.{i+i)h f. 

J2 du [X{u) - X{ih)]dz 

i=0 ■'ih J 

i-1 Ai+l)h n 

/ {u — ih)du i X'{ih)dz 
i=o -^^^ 

+ J2 . {u-ihfdu {I - 5) I M' {s) \s^s, dzdddu (131) 

„-_n J ih Jo J 



where X{u) = p{0,u, x, z) ^ Xv{w, v)d'w^ D^H{u,jh, z, v), Si = ih + S{u — ih). As in the treatment of 
Ti, we obtain that 



/ {u — ih)du I X!{ih)dz 

h "'^^ /'•''* f 
= -^/i / Xu{s,v)ds i p{0,ih,x,z)D''^Ao{ih,jh,z,v)dz 

^ i_n Jih J 



=0 

oXl^ / / p{0,ih,x,z)D'^Hi{ih,jh,z,v)dz 

^ i=o -^'^ 

2 y p{0,ih,x,z)D''^H{ih,jh,z,v)dz, (132) 



where 



Ao(s, j/i,i;,D) = (i -2LL + L )p{s,jh,z,v), 



I ^ ( d'-(Tij{s,z) d^aij{s,v)\ d'^p{s,t,z,v) 
m{s,t,z,v) = (^^-^ ^-^J 



dzidzj 



for I = 0,1,2 with iJo = H. The differential operator Aq was introduced before equation (47). It is a 
fourth order differential operator. From the structure of this operator and from (132) it is clear that it is 
enough to estimate 

,71-1 ^jh , 

^-oY^^ Y.^. xAs,v)ds / p{0,ih,x,z)D''+^'p{ih,jh,z,v)dz^jh,T,v,y)dv (133) 
j=0 •' i=0 •^^'^ 

for = 3, = 3. To estimate (133) we consider three possible cases: a) jh > T/2,ih < jh/2 =J> 
jh -ih> T/4 b) jh > T/2, ih > jh/2 =^ ih > T/4 c) jh <T/2=>T - jh > T/2. In the case a) we 



48 



apply integration by parts. This transfers two derivatives to p{0,ih,x, z). This gives 

SE'^ r xAs,v)ds [ Dl''+-'p{0,ih,x,z)D-',+^'-^''-^'p{ih,jh,z,v)dz 
Jih J 



Jo u^-^ {jh-uy-^^^"- ' 

<C{e)h'-^%3hr-'4>^{v-x) 



and 



|/| < C{e)h 



l-2e [ 



du 



< C(£)/i3/4TV2n-(V4-2e-V2)^^(y _ a;) 

< C{s)Ty^-'h^/^+'<j>^iy-x), 



(134) 



where (5 = (1/4 - 2e - >i/2) > if x < 1/2 - 4£, < e < 0, 05 (see the condition (B2)). In the case b) 
we apply integration by parts and transfer four derivatives to p{0,ih,x,z). This gives the same estimate 
as in (134). At last, in the case c) we make an integration by parts and transfer three derivatives to 
^{jh,T,v,y) and one derivative to p{0,ih,x, z). This gives the same estimate as in (134). To pass from 
D^p{ih,jh,z,v) to Di^p{ih, jh, z.v) we use the following estimate 

\D''J{ih, jh, z, v) + Dlip{ih,jh, z, v)\ < C<l)^^jj^{v - z). 

Clearly, the same estimate (134) holds true for the other summands in the right hand side of (132). This 
gives 



o "12^ ll2^ X,.(s,i^)c^s / piO,ih,x,z)D'^Hi{ih,jh,z,v)dz^{jh,T,v,y)dv 

<C{e)T^'^-'h^'^+'ct>^{y~x), 

2 



/ xAih,v)^h / p{0,ih,x, z)D^H{ih, jh, z,v)dz^{jh,T,v,y)dv 

i=0 '' i=0 '' 



<C{s)T'/'-'h'/^+'^^{y-x). 
We now estimate the second summand in the right hand side of (131). Similarly as in (50) we obtain 



dzdSdu 



^ / {u- ihfdu I (1 - (5) / X"{s) I 

^-QJih Jo J 

= X1 /. {u-ihf / (1-5)^ / xAT,v)dT / p{0,s,x,z)D^Ak{s,jh,z,v) |s=s, dzdSdu 
~X] / {u-ihf I {l-5)xAs,v) I p{0,s,x,z)D^Ao{s,jh,z,v) |s=s, dzd6du 

i—Q-Jih Jo J 

J-1 r(i+l)h pi p 

/ {u-ihf I {l-d)x^{s,v) p{0,s,x,z)D^Hi{s,jh,z,v)\s=si dzdSdu 

i=0 "^0 

-V/ {u-ihf {l-S) ' p{0,s,x,z)D''^H{s,jh,z,v)\s=s, dzdSdu, (135) 

.■n Jih Jo J 
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where the operators Ai,i = 1,2,3,4, are defined as follows: 

Ai{s,jh,z,v) = {L^ -31^1 + 311"^ -L^)p{s,jh,z,v), 

A2{s,jh,z,v) = {LiH + 2LHi){s,jh,z,v). 

A3{.s,jh,z,v) = [{L-L)Li+2{Li-Li)L]p{s,jh,z,v), 

A4{s,jh,z,v) = H2{s,jh,z,v). 

The operator Ai was introduced in (51). For this operator it is enough to estimate for fixed p, q, r, I 

As in (56) we obtain that this term does not exceed 

C{e)h^/'-%jhf^-'cl>^{v-x). (136) 

It follows from the explicit form of these operators that the same estimate (136) holds for ^2, ^3 and A4. 
The other three terms in the right hand side of (135) do not contain the factor j'''^'' Xu{'''i'")d'r and they 
can be estimated separately. Clearly, it is enough to estimate the term containing Aq. The remaining 
two summands are less singular. From the explicit form of Aq (compare also (46)) we obtain that it is 
enough to estimate for fixed q,l,r 



(u-thf I {1-S)x.{s,v) I p{Q,s,x,z)Dl (^^^^^^Ih^ ) {s,jh,z,v) dzdSdu. 



Analogously to (47) we get that this term does not exceed 

C{s)h'-'%jhr-'ct>^{v-x). (137) 
Now from (131), (134), (135), (136) and (137) we obtain that 

\p(»h{L- L)ni] <»h $(0, T,x,y)-p® [{L - L)ni (»'^ $](0, T, x, y) 

< Ch^^^+',p^{y - x) (138) 

for some ^ > 0. The other replacements can be shown analogously. Thus we come to the following 
representation 

Ph{0,T,x,y) -p{0,T,x,y) 

= Vh [ni $(0, T,x,y)+p® (SRi $)(0, T, x, y)] 

+h [n2 $(0, T,x,y)+p(3 (K2 K *)(0, T, x,y)+p (3h (5?3 K *)(0, T, x, y)] 
+h [ni $ + p (5?i $)] 0'^ (5?i 0'^ $)(0, T, X, y) 

+ {Ll - L^)p{0, T, X, y)-^p® {L' - L')p{0, T, x, y) + 0{h'+'C^{y - x)). (139) 

We now further simplify our expansion of ph — P- We start by showing the following expansion 

Ph{0,T,x,y) - p{0,T,x,y) 

= Vh{p .Fi[pa])(0, T, x,y)+h{p(^ ^2bA]) (0, T, x, y) 
+/i (p J^i b ^1 [pa]]) (0, T, X, y) 

+ ^P {Ll - L^)p{0, T, X, y)-\p® [L' - L')p{0, T, x, y) + 0{h^+'C^{y - x)), (140) 
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where for s G [0,t - h],t G {h, 2h, T} 

PA{s,t,z,y) = {piSi'h^){s,t,z,y) 

= Pis,t,z,y)+ ^ h j p{s,jh,z,v)^i{jh,t,v,y)dv. 

s<jh<t—h 

Here = H + H igi'h H + H iSi'^ H iSi'f^ H + .... We now treat the term p (g) Z^i(s, t, x, y). 
p'SiL7ri{s,t,x,y) = f dr f p{s,T,x,v){t - t) ^ ^"^'^'^'^^ Dl(Lyp{T,t,v,y))dv 

■'^ W\=3 



= -^-]fdv [ p{s,T,x,v)(f Xu{u,y)du\ ■^D''^p{T,t,v,y)di 



kl=3 

= / + //. 



kl=3 



(141) 



By integrating by parts w.r.t. the time variable we obtain for /. 

I = -T.^ [dv \p{s,T,x,v) ( fx,{u,y)du) D^J{T,t,v,y) 

- DyP{T,t,v,y) ^ ^^(^'^^^'""^ J xAu,y)du-p{s,T,x,v)xA'^>y)^ dr 
= -] I dv p{s,^^,x,v) I / xA'U;y)du \ D^p{^-^,t,v,y) 

ii/i=3 ^' L ) 

+ m 7i f / Xv{u-,V)du \ D';;p{s,t,x,y) 



1 f 2 



kl=3 



+ ^ — / c^'T / X,.(w,?/)du) / L^p{s,T,x,v)Dlp{T,t,v,y)dv 



~ X] ~!_^ xAT^y)dT j p{s,T,x,v)D^p{T,t,v,y)dv. 



kl=3 

For the second term we get 



(142) 



+ ^ -J (drf I x>.(w,y)(iuj / L^p{s,T,x,v)D''^p{T,t,v,y)dv 



(143) 
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From (141)- (143) we have 

piSiLTTi{s,t,x,y) = ni{s,t,x,y) -|-p(g) LTTi{s,t,x,y) - pigi Mi{s,t,x,y). 

This shows that 

Tri{s,t,x,y) +p(g)'Sii{s,t,x,y) = TTi{s,t,x,y) 

-\rp (g) Ltti {s,t,x,y) - p® Ltti {s,t,x,y) + p ® Mi{s,t,x,y) - p ® Mi{s,t,x,y) 
= p0Mr{s,t,x,y). (144) 

It follows from (144) and the definitions of the operations (8) and (8)^ that 

Vh [ni 4>(s, t, X, y) + (p (g) 3?i) O'^ «>(s, i, x, y)] 
= Vhiwi + p (8) 3?i) (E)^ $(s, i, a;, y) 
= Vhip O Ml) $(.s, t, X, y) 

= \/h ^ h I {p® Mi){s,jh,x,z)^{jh,t,z,y)dz 



0<jh<t-h 



^{jh,t,z,y)dz 

^{jh,t,z,y)dz 



= Vh h / du p{s,u,x,v)Mi{u,jh,v,z)dv 

0<jh<t-h J [Js J 

= Vh h / dux[s,jh] / p{s,u,x,v)Mi{u,jh,v,z)dv 

0<jh<t-h J iJs J 

= Vh f du f p{s,u,x,v) ^ •^"^J,'^^ X Dj; ^ /ix[s,i^] / p{u,jh,v,z)^{jh,t,z,y)d 

•'^ \v\=Z ^' 0<jh<t-h ■' 

du J p{s,u,x,v) X ^ — D'^pA{u,t,v,y)dv 



dv 



IH=3 

= \/^(p8)jri)[pA](s,i,x,y). (145) 

Here, xi^^jh] denotes the indicator of the interval [s,jh]. Using similar arguments as in the proof of (145) 
one can show that 

h [n2 ®'h 4'(s, t, X, y) + {p(E> ^2) <E>'h ^{s, t, x, y) +p (E>h (3?3 K t, x, y)] 

= h{p®J^2)\pA]{s,t,x,y) + hp®Ui ®'^^{s,t,x,y). (146) 

For the first two terms in the right hand side of (139) we obtain from (145) and (146) 

Vh [ni $(0, T, X, y) + {p®^i) $(0, T, x, y)] 

+h [^2 ®'h $(0, T,x,y)+p® (3*2 ®'u $)(0, T, x,y)+p ®h (3?3 ®'h 3')(0, T, x, y)] 
= Vh{p J^i)\pa]{0, T, x, y) + h{p T2)Ipa]{s, t, x, y) + hp® Hi (g)'^ $(s, t, x, y). (147) 

Using (145) we get 

/i 1 ® $ + p ® (3?i $)] ^U^R 1 ® L (0, T, X, y) 
= h{p ® Ti[pa]) ®'h (3?i ®'h *)(0, T, X, y) 
= hp®Tt [pa ®'h (3*1 < *)] (0, T, X, 2/). 
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Note that 



hp Hi (g)'^ $(s, t,x,y) = h j^du j p{s, u, x, v) ^ '^"^^^'^^ D'^[wi y) 



= /ip0^i[7ri 0'^ $](s,i,a;,y). 



For the proof of (140) it remains to show that 



hp J^ilm ®'t, $ + PA ®'h {^i ®'h *)](0, T, X, y) 

= h{p®T^\p® Ti\pa]\) (0, T, X, y) + 0{h^+\^{y - x)). (148) 



We will show that 



hp (g> ^1 [(p - pa) ®'h $)] (0, T, a;, y) = 0(/ii+^Cv^(y " a;)), (149) 

hp®Ti\p (3?i $)] (0, T, y) - /ip ^1 [p (3?i $)] (0, T, a;, y) 

= o(/ii+*CvT(y-^))- (150) 

Claim (148) follows from (149), (150) and (145). The estimate (150) can be shown similarly as in the 
proof of (138). An additional singularity arising from the derivatives in the operator can be treated 
by using the additional factor h in (150). To estimate (149) note that from the definition of and $ 

|(5ii0'^$)(j/i,T,a;,y)| <C(£)/i-^(T-i/i)-VvTr77r- (151) 

Then we use the fohowing estimate which can be proved by the same method as in the treatment of Ti , 
where an estimate for {p — p'^){Q, jh,x,y) was obtained. 

\{p-PA){u,jh,v,z)\<Ch^'^4>,^jj^{z-v). (152) 

From (151) and (152) 

\{p - pa) (^'h (3?i 0k *)(w, T, v,y)\< C{e)h^'^-'{T ^ uf^p^r^iv - v) (153) 

For an estimate (149) it is enough to estimate a typical summand of the sum of the detailed representation 
of the left hand side of (149). E.g., for \v\ = 3 we have to estimate 

h I du I p{Q,u,x,v) ^"^^^^^^ D^^l ^ h I {p - pA){u,jh,v, z) 

° {j-u<jh<T-h} ■' 

x(5Ri ®'h^){3h,T,z,y)dz\dv 

fT/2 



= h ... + h . 

Jo Jt/2 

= 1 + 11. 



For an estimate of II we apply integration by parts and transfer three derivatives to p{0,u,x,v) "^^f' 
Using (153) we obtain the fohowing estimate 



\II\ < C{e)h^/^-' I ^^L^du(t>^{y-x) 
Jt/1 ' 

< C{e)h^'^-'T'cl)^{y-x). (154) 
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For the treatment of I we consider two cases: a) jh — u > T/4 and b) jh — u < T/4 T — jh > T/4. 
Similarly as in (42) the difference h{p — pa) can be represented as 

h{p - PA){u,jh,v,z) = h{p iSi H - p ig)'^ H){u,jh,v,z) 
+h{p ®H ~p®'j^H) jh, V, z) 

u, T, V, z')H{t, jh, z', z)dz' 

/ dr I {\{T,z')-\{ih,z'))dz' + h{p®H -p®'f^H)®'f^^i{u,jh,v,z) 
Jih J 



(155) 



where X{t,z') =p{u,T,v,z')H{T,jh,z',z),^\{ih,jh,z'z) = H{ih,jh,z'z) + H®'f^H{ih,jh,z'z) + ...,j* = 
j*{u) = [^] + 1. Here [x] is equal to the integral part for noninteger x and equal to a; — 1 for integer x. 
For /', case a), we have jh — t > T/5 for n large enough. With the substitution v + v' = z' we obtain 

D'^h / dr p{u,T,v, z')H{T,jh, z' , z)dz' 

D'^h J dr j p{u,T,v,v + v')H{T,jh,v + v' , z)dv' 



< Ch 



dr 



u {jh - t) 



<Cn-'^^jj^{z-v)=CT'h'cP^jj^{z-v) (156) 
For the proof of (156) we used the following estimate from Friedman (1964) (Theorem 7, page 260) 



\D';;p{u,T,v,v + v')\ < C{t - u)-"^^^ exp 



c\ 



— U) J ' 

For J I'{u,jh,v, z){^i ®'i^ ^){jh,T, z,y)dz, case b), it is enough to estimate for \u\ = 3 

^,d^piT,jh,V + v',z)_^^, 



hi dr \p{u,T,v,v + v'){aik{T,v + v') - aik{T,z))]D'^,- 
x(SRi ®'f,^){jh,T,z,y)dz. 



(157) 



For an estimate of this term we transfer five derivatives from p to (3?i ®'i^ ^){jh, T, z, y) and we use the 
following estimate for = 5 

\D'i,p{T,jh, V + v', z) + D^p{T,jh, V + v', z)\ 
<C{jh-T)-''/^^jj^{z-v-v'). 

We obtain that (157) does not exceed 



C{j*h - T)hT 



-7/2, 



^{y-v) < Ch'+'T'-'n-('-'-'-/^UvT=^{y-v) 



o{h'+'T'-' 



(158) 



We used that for any < 5 < 1 it holds that k < ^ , see condition (B2). For an estimate of 
/ IV (u, jh,v , z){^\ ®'y^ ^){jh,T, z,y)dz we use the decomposition (43). For getting an estimate for the 
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terms in //' that contain the first derivatives X'{ih, z') we use the identity (45) and similar arguments as 
already used in the estimation of / 1' {u, jh,v, z){^i (8)'^ <^){jh,T,z,y)dz. The estimate for terms in //' 

containing second derivatives A {ih,z') follows from (53) and (54). Finally, for III' the same estimates 
hold because of smoothing properties of the convolution ...(S)^$i(«, j/i, v, z). This implies (149) and, hence, 
the expansion (140). 

Asymptotic replacement of pa by p. Now, we compare hp^ J^2\pA]{0,T,x,y) with hp J^2\p]{0,T, x,y) . 



Note that for 2>!: < (5 < |, = 4 

h j du j p{0,u, X, z)x^{u, z)D'^p{u,T, z,y)dz 



<Ch^+'{T-h')-^<P^{y-x) 
< C7/i1+(^-2-)t2^<^^(j/ - a;), (159) 
< C/i^+(*-^")r2^(^^(t/ - x). (160) 



(Tn^ - n^/i'5)2 
h du D''^\p{0,u,x, z)x^{u, z)]p{u,T, z,y)d 

JT-h^ J 

The same estimates hold for pa{u, T, z, y). Hence, it sufltices to consider u G [/i*, T — h^]. We now treat 

h j^^ du j p{0,u,x,z)Xt,{u,z)D'^{p-pA){u,T,z,y)dz 



l-T/2 pT-h" 
= h ... + h 

Jh^ Jt/2 



...=I + II. 



(161) 



By using (152) we get 



h f du f D",\p{0,u,x, z)xAu, z)](p - PA)iu,T, z,y)d 

Jt/2 j 

< Ch^'^n''(t)^{y - x) = Ch^^^-^T^cp^iy - x) 
= Ch'+'^cP^{y~x),^>0 

For u e [h^,T/2] it holds that 

T/2 



(162) 



< 



Ch^/^-'<j>^{y-x). 



Note that the condition <5 < | implies that 3/2 — ^ > 1. It follows from (159)-(163) that 



hp ® J^sbA] (0, T, X, y)-hp® J-2[p](0, T, x, y) = 0{h^+'' <j)^{y - x)). 



(163) 
(164) 



For the proof of 



hp(g)J'i[p(S>J^i[p-pA]]=0{h^+^<j>^{y-x)) 
we consider a typical summand for fixed ly, {vl = 3, 



du / p{0,u,x,z)xAu,z)D^z 



j j p{u,T,z,v)xAT,v)Dl{p-pA){T,T,v,y)dv 



dz. (165) 
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As before it suffices to consider the integrals for u G [h^, T — h^]. The integral in (165) is a sum of four 
integrals 



h 
h 

h 



h du ...D^ / dr 

Jh^ J Ju J 

/■T/2 f r 
h du ...£>^ / dr 

Jh^ J J(T+u)/2 J 



)/2 

pT-h' r(T+u)/2 I- 

= h du ...Dl / dr / 

Jt/2 J Ju J 



T/2 

-T-h 

T/2 



du / ...D': 



dr 



(T+u)/2 



(166) 



Note that in the integrand in I2 it holds that t — u> T/4. By applying integration by parts w.r.t. v and 
(152) we get 

I/2I <C/i3/'-"T>^(2/-a;)). (167) 

Furthermore, in the integrand in /4 it holds that u > T/2, t — u> /2, T — u > . Using integration 
by parts w.r.t. z we obtain 

I/4I < C/l3/2-*T-V2<^^(y _ :r)) < CT-/2/i3/2-W2-5^^(y _ ^))^ (168) 

where, by our choice of 6, 3/2 — k/2 — S > 1. For an estimate of I3 we use the representation 

(p-PA)(r,T,v,y) = {p®H~p®'f^H){T,T,v,y) 
+{p®H-p ®; H) $i(r, T, V, y) 

= j ds j p{T,s,v,w)H{s,T,w,y)dw 

+^[p<8>'AHi+Ao)](T,T,v,y) 

n-l /.{i+l)h 



+ H / {t-ihf j {l--f)y^ I p{T,s,v,w)Ak{s,T,w,y)\s=ih+'y{t-ih)dwd^dt, 
^ i=j* -l^h Jo J 



+ {p(»H-p (g)'^ H) (g){, $i(t, T, V, y), 



(169) 



where j* = j*{T) = [r//i] + 1. As above [x\ denotes the integer part for nonintegers x and it is equal to 
X — 1 for integers x. The quantities Hi and ^fe, fc = 0, 1, 2, 3, 4, have been defined (52) and 

^i{ih,ih,z,z) = H{ih,ih,z,z') + H (»'f^ H{ih,i' h, z, z ) + .... 

To estimate -D^(p — Pa){t, T, v, y) we note that 

DJ^ / ds I p{T,s,v,w)H{s,T,w,y)dw 



= h 



< Ch 



D'^ j ds I p{t, s,v,v + w')H{s,T,v + w' ,y)dw' 
ds 



{T-sy 



2-25, 



{y - v). 



(170) 
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Furthermore, 



D:ip(g>'^H,]{T,T,v,y) 



T<jh<T-h 



< 

- 2 



J2 h I D'^\p{T,jh,v,v + w')Hi{jh,T,v + w',y)]dw' 



'<jh<T-h^ /2 



^ E E ^ [ DT'^"' W, jh, v,v + w'Mjh, T,v + w', y)dw' 

i,k=lT-h^ /2<jh<T-h 

< Ch^-^^^^^iy -v) + Ch^-^^^^(f>^^{y - v). 
Because of the structure of the operator Aq it is enough to estimate for fixed i, I, k 

,^d^p{jh,T,v + w',y)^ , 



dw'-dw', 



-dw'. 



T<jh<T-h 

With the same decomposition as in (171) we obtain that (172 ) does not exceed 

Ch'-'^/'<P^{y-v). 

From (171) and (173) we obtain that 

^ mp^'f, iH,+Ao)]iT,T,v,y)\ < Ch^+'^^^iy-v) 



(171) 
(172) 

(173) 
(174) 



for some 7 > 0. It remains to estimate the last summand in (169). It follows from the structure of the 
operators Ak, fc = 1, 2, 3, 4, that it is enough to estimate 



X / {t-ihf {i-j)j2 d: 



p(r, s,v,v -\-w ) 



dw^dw'idw'pdiv'q 



for fixed i,j,p,q. As above, we obtain that (175) does not exceed 

„1 /.I n-l 



ih+i{t-ih) dw'd-ydt 
(175) 



r r 1 
Jo Jo [{^h-T}- 

= Ch^<f>^{y - v) f j\l - 7) 

""^ " {i:3*h<ih<T+h^/i} 

+Ch'<t>VT^{y - v) f z^ /\i - 7) E 

= 1" + 11". 



+ 7/12;]^/^ [(n — ^z)h — ih\ 



-d'ydz 



(176) 



(177) 



Now, 



d'ydz 



{i:j*h<ih<T+h^/4} 



< Ch^-^h-^'/^./j^iy-v) f\^-' [' ^-^—^ y h— ±—-—djdz 

'Jo Jo r [{ih-r)+'yhz]^-- 



{i:j*h<ih<T+h^ /4} 



< C{s)h'-'-''/^cP^{y-v). 



(178) 
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Using inequality {h — 'yz)h — ih = {n — i)h — 'fzh > h{l — 72:) > h{l — 7) we obtain that 
\II"\ < Chh-''"^^{y - v) f d-i ^ 

"'^ "'^ {i:T+h^/i<ih<T-h} 

< Ch'-''/^<l>^{y-v). (179) 
Now from (169), (170), (174), (175), (178) and (178) we obtain that 

\D:ip(^ H - p(g>'^ H){T,T,v,y)\ < Ch-'^^j^iy - v) (180) 

for some positive 7. The last summand in the right hand side of (169) admits the same estimate (180) 
because of the smoothing properties of the operation (Xijj. Hence, 

\D-^{p-pA){r,T,v,y)\<Ch''^^{y-v). (181) 

Making the change of variables v = z + v' into (165) we get that the integral w.r.t. v is equal to 



j dr j p{u,t,z,z + v')x„{t,v)D'^{p~ pA){T,T,z + v',y)dv' 



(182) 



Taking into account (181) and applying integration by parts in (182) we obtain that (182) does not exceed 
From (165) and (183) we obtain that 

\h\<Ch'+^ct>^{y-x). 
The estimate for I\ can be proved analogously to the estimate for 73. Thus, we proved that 

hp®Ti\p®Ti\p-pA]] = 0{h}+^4>^{y-x)). 

The estimate 

h^'^p ® ^1 [p - pa] = 0{h}+^(l>^{y - x)) 
can be proved by using the same decomposition of p — pa- This completes the proof of Theorem 1. 
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